ON IMPACT DAMPERS FOR NON - LINEAR AND 
SELF - EXCITED VIBRATING SYSTEMS 


by 

SHYAMAL CHATTERJEE 


ro p* 

o 

CHF> 

ON 



DEPARTMENT OF MECHANICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

August, 1995 



On Impact Dampers for Non-Linear and 
Self-Excited Vibrating Systems 


A Thesis submitted 

in partial fulfilment of the requirements 
for the degree of 

Doctor of Philosophy 


by 

Shyamal Chatterjee 


to the 

DEPARTMENT OF MECHANICAL ENGINEERING 
INDIAN INSTITUTE OF TECHNOLOGY 
KANPUR 
August 1995 


Dedicated ta 


(TUj^ &el(Med 


ponenixx and ^uAten 



CERTIFICATE 


This is to certify that the present work titled ON 
IMPACT DAMPERS FOR NON-LINEAR AND SELF-EXCITED VIBRATING 
SYSTEMS by Mr. Shyamal Chatterjee has been carried out under 
our supervision and that, it has not been submitted elsewhere 
for a degree. 


CUjul 

(A. Ghosh) 



(A. K. Mallik) 


Professor 


Professor 


Department of Mechanical Engineering 
Indian Institute of Technology, Kanpur 



ACKNOWLEDGEMENTS 


It is a rare piece of fortune to have Professor A.K. Mallik 
and Professor A. Ghosh as my thesis supervisors. Only 'thank you' 
would be too cold a word to express my deep sense of gratitude. So 
let it be never expressed but felt for ever. 

My sincere thanks are due to Professor N. Vays of our 
vibration lab. 

I specially thank Professor A.K. Mallik for introducing me to 
the field of Non-linear Vibrations. Frequent visits of Professor 
J.K. Bhattacharyya to our laboratory were a continuous source of 
inspiration for our studies on CHAOS. I would also like to thank 
Professor M.R.M. Rao for his advanced courses on O.D.E. and 
Liapunov stability. I am very happy having Professpr V. Gupta as 
one of my teachers. His friendly attitude to us can never be 
forgotten . 

What, if I didn't have Ravindra, Rajiv, Chandrashekhar around 
me! It would have been simply impossible for me to live and work 
in the cool of our vibration lab. They kept it warm. They all were 
the sources of vibration, energy, fun, enjoyment and so on. 

Soumya, Anirvan, Biswajit, Indra, Sudipta, Rahul, Biswaroop, 
Atanu and others brought chaotic order in my hall four days. It 
was they who made my six years stay at IITK, so short. 

My special thanks are due to Mr. M. M. Singh for his kind 
assistance. 

And my little, window-side table - you never felt tired 
bearing my load and sharing my frustrations, my pains, my joys. 
How do I acknowledge you! Or how do 1 express my love to you, my 


V 


cute, red cycle for carrying my weariness along the asphalt roads 
of IITK. I enjoyed your company on the way of my journey from the 
glass-doors of the C.C. to the bed of my dreams in so many dark 
silent nights. 

Thanks every one I know. 


y/u/ama£ 



SYNOPSIS 


Linear/non-linear vibration neutralisers (absorbers) , 
designed for controlling the vibrations of lightly damped systems 
are available in various forms. The hardware design of vibration 
neutralisers consists of a secondary mass connected to the primary 
vibrating body by a linear elastic member. It is known that a 
non-linear coupling force between the primary and secondary masses 
greatly enhances the effectiveness of a vibration absorber. An 
impact damper, in the form of a loose secondary mass provides a 
clearance type non-linear coupling force. This clearance 
eventually gives rise to repeated impacts between the primary and 
secondary masses and the vibration is attenuated mainly through 
the transfer of momentum. The simple constructional features of an 
impact damper make it very attractive for use in real life 
systems. Impact dampers have been successfully used to control the 
vibrations of various linear systems with well separated modal 
frequencies. The theoretical design methodologies of an impact 
damper are also developed based on the assumption that the primary 
system is linear. However, the non-linearity present in the 
primary system tends to show up its prominence when the system is 
strongly excited and needs to be controlled. Moreover, the 
oscillatory characteristics of certain systems, such as 
self-excited systems, are fundamentally non-linear. The major 
objective of the present thesis is to study the performance and 
design of an impact damper for non-linear systems. Three different 
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types of models of non-linear systems viz., 

(i) harmonically driven Duffing' s oscillator, 

(ii) autonomous van der Pol oscillator and 

(iii) harmonically excited van der Pol oscillator 

are considered in details. Besides the design of an impact damper, 
studies on various complicated dynamical responses of the systems 
under consideration are also carried out in the light of the 
bifurcation theory and chaotic dynamics. A brief outline of the 
contents of the thesis is given below. 

The systems under consideration can be classified as 
piecewise smooth, non-linear systems. In chapter 2, the equivalent 
linearization procedure has been extended for such systems and is 
termed as the equivalent piecewise linearization (EPL) method. The 
efficacy ,of the method, as compared to the simple harmonic balance 
method (HBM) , for piecewise smooth, non-linear systems has been 
demonstrated by considering a Puffing's oscillator with motion 
limiting stops (Vibroimpacting oscillator) . The stability analysis 
of the solutions is also presented. All the analytical results are 
verified by numerically integrating the equations of motion. 

In chapter 3, the design and performance of an impact damper 
for controlling the high-amplitude vibration of a harmonically 
excited. Duffing type primary system are discussed. Again, the 
results are obtained by the HBM, EPL and numerical integration. A 
detailed parametric study is carried out to demonstrate the 
effects of various parameters. The procedure for arriving at an 
optimal design is presented. The question of controlling the 
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subharmonic response has also been addressed to. Finally, the 
possibility of using the impact damper discontinuously in an 
active manner has been demonstrated. In this mode, the impact 
damper is shown to act as a migration controller which changes the 
initial conditions so as to bring the original system within the 
basin of attraction of a lower amplitude solution. Thus the 
feature of multiple stable solutions of a non-linear primary 
system is put to use. 

In chapter 4, the performances of an impact damper towards 
controlling the autonomous (limit cycle amplitude) and 
harmonically entrained forced oscillations of self-excited 
oscillators, are discussed separately. In the latter case, after a 
detailed parametric study, an empirical relation for the optimum 
design is established. A sensitivity analysis of the design with 
reference to the self-excitation and forcing parameters is 
included.' 

In chapter 5, a study of bifurcations and chaos in 
self-excited, autonomous oscillators with an impact damper is 
presented. Results have been obtained by considering various 
models of the self-excitation mechanism. The major emphasis is on 
the detection of the types of bifurcations which are independent 
of the exact nature of the model. In this chapter, a numerical 
technique for analysing the eigen properties of the numerically 
integrated response has been proposed. 

In chapter 6, different routes to chaos for a forced, van der 
Pol oscillator with an impact damper have been investigated. Some 



uncommon intermittency routes to chaos have been detected. Finally 
in chapter 7, the overall conclusions and some suggestions for 
future works are included. 

The major conclusions of the present thesis are listed below. 

(i) The proposed EPL technique contains information in both time 
and frequency domains and consequently gives better results than 
the HBM which is simply a frequency domain technique. This method 
is particularly suitable for obtaining the time domain response of 
a class of piecewise smooth, non-linear systems. 

(ii) An impact damper is effective, both in passive and active 
on-off forms, for controlling the vibration of non-linear systems 
having Duffing type restoring force characteristics. Both the 
primary and secondary resonances can be controlled. A simple 
formula is obtained for providing an initial design and it has 
been shown how to improve upon to obtain a final design, 

(iii) An impact damper can be used to control both free and forced 
self-excited vibrations. Sensitivity of the designed performance 
towards the inaccuracy of the estimated model parameters is to be 
incorporated in the design rules. An empirical design formula for 
an impact damper attached to a forced van der Pol oscillator is 
provided after a detailed parametric study. The role of the 
coefficient of restitution for various situations has been clearly 
brought out. 

(iv) Various instabilities occur in an autonomous self-excited 
system with an impact damper. The gross bifurcation superstructure 
of the systems under consideration is more or less independent of 
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the exact nature of the model of the self-excitation. A numerical 
technique has been developed for analysing the eigen properties 
of the stable solutions. 

(v) The dynamics of a harmonically driven van der Pol oscillator 
with an impact damper is rich with different types of 
intermittency catastrophes viz., type II, type III etc. A new type 
of intermittency after a subcritical, symmetry-breaking pitchfork 
bifurcation is seen to exist in the system under consideration. 
Several instabilities of the system has a close resemblance to the 
hydrodynamic instabilities such as Rayleigh-Benard convection. 
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Figure 5.20 Bifurcation diagram of the periodic solutions 

with two symmetric impacts per cycle for the form 
of the function f shown in Figure 5.2(b) with 
various values of r^. x = 0.4, e = 0.25. a) o two 
impacts per cycle and • four impacts per cycle for 
r^= 0.05; A two impacts per cycle and A four 

impacts per cycle for r^ = 0.075; o quasiperiodic 
motion; b) o two impacts per cycle and • four 

impcts per cycle for r^ = 0.08; ® two impacts per 
cycle for r^ == 0.085; o quasiperiodic motion; c) o 
two impcts per cycle and • four impacts per cycle 
for = 0.09; d) • four impacts per cycle , no 
symmetric two impacts per cycle exists for r^ = 

0.18. IW represents interwindow space. 

Figure 6.1 Bifurcation set in r^- d^ plane, n = 1, X' = 0.7, 

Fq= 0.5, c = 0.1. The zone enclosed between the 

lines AB and CD corresponds to stable symmetric 
motion with two impacts per cycle. AB and CD 
correspond to supercritical and subcritical Hopf 
bifurcations, respectively. 

Figure 6.2 Bifurcation sets in e - d^ plane, n = 1. a) Fq= 

0.5, r^= 0.25; b) Fq= 0,5, = 0.7; c) = 0.7, 

r = 0.25. Hopf bifurcation (the upper line 

corresponds to supercritical bifurcation and the 
lower line corresponds to subcritical 
bifurcation) , - - - pitchfork bifurcation 

(subcritical) . 

Figure 6.3 Bifurcation set in F^- n plane, x = 0.7, c = 0.1. 

AB - supercritical Hopf bifurcation, BC 
subcritical Hopf bifurcation, CA - subcritical 
pitchfork bifurcation : d^^ 1,6, rj^= 0.45. EF - 

supercritical Hopf bifurcation, EGF - subcritical 
Hopf bifurcation : d^ = 2.9, r^^* 0.25. 
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CHAPTER-1 

INTRODUCTION 


1. 1 Introduction 

A rich variety of qualitatively different phenomena may be 
exhibited by a non-linear system. Such complex behaviour can never 
be predicted through the corresponding linearized system. In the 
last fifteen years, there is an explosion of literature related to 
the study of non-linear systems in various branches of physics and 
engineering. While the mathematician and physicists are attempting 
to analyse and categorise various non-linear phenomena, the 
engineers are trying to use these phenomena towards improving the 
performance of practical systems. 

A mechanical vibrating system is often modelled through a 
combination of spring, mass and dashpot where the spring 
represents the restoring force and the dashpot represents the 
dissipative force. If the restoring and/or the dissipative forces 
present in the system are represented by some non-linear functions 
of the motion variables, then the system is called non-linear. 
Mostly these functions are considered to be smooth in nature. Such 
smooth non-linearity may arise due to large motion, material 
non-linearity and coupling of different modes of motion etc. But 
there are situations where an abrupt change in the restoring 
and/or dissipative forces or in the state of the system takes 
place. Such an abrupt change automatically brings in non-linearity 
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in the system, even if the nature of the restoring and dissipative 
forces basically remain linear. Speaking in mathematical terms, 
such non-linearity is introduced through the non-analytic nature 
of the restoring and dissipative forces which appear in the 
differential equations describing the model of the system. Such a 
system is often called a non-dif f erentiable or piecewise smooth 
system. 

The latter kind of non-linearity is present in a variety of 
man-made devices which engineers cannot afford to overlook. 
Connected components with unavoidable clearances in the joints, 
machining conditions during intermittent cutting, backlash in 
meshing gears, rotor-bearing mounting with clearance, colliding 
machine components such as impact hammers etc., are some common 
examples in this category. Sometimes clearances and impacts are 
deliberately introduced for improving the system performance. In 
the field of vibration control, the typical examples included the 
motion limiting stops such as elastic fenders for moored ships, 
vibration arresters for heat-exchanger tubes. 

The use of a secondary system, elastically connected to a 
primary system, for controlling the near resonance, high amplitude 
oscillation of the latter is quite common. Such a secondary system 
is normally referred to as a dynamic vibration neutraliser (DVN) . 
The analysis and design of such DVN's based on linear theory, in 
various forms, are widely reported in the literature on 
engineering vibrations. The idea of using a non-linear DVN for 
better effectiveness has also been proposed. But in practice, the 
introduction of a desired and controllable non-linearity in a 
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system is by no means an easy task. It is possibly the clearance 
which can generate non-linearity in a system in an easy and 
controllable fashion. Based on this idea, an impact damper is 
conceived where a loose secondary mass undergoes repeated 
collisions with the primary system and in the process reduces the 
unwanted vibration of the primary system. 

Though the non-linearities may be effectively used for the 
purpose of vibration control, one must also remember that a 
non-linear system may become unstable in a number of ways to mar 
the system characteristics. This aspect of stability deserves 
serious attention at least in a qualitative sense. Recent studies 
of qualitative non-linear dynamics have shown ample promise to 
throw light on a number of complicated behaviours frequently 
observed in laboratory experiments and practical situations. For 
example, the noise in an apparently simple looking system (working 
in a more-or-less deterministic dynamic environment) may have its 
origin in the simple non-linearity provided by the clearance. 

One of the major objectives of the present thesis is to 
investigate the role of an impact damper in controlling the 
vibration of inherently non-linear systems. Design philosophies 
and methodologies have been discussed for a class of non-linear 
oscillators such as Duffing' s oscillator and van der Pol 
oscillator. Besides the aspects relating to vibration control, the 
characteristics of the complex behaviour exhibited by such 
non-linear systems with clearances and impacts have been studied 
in the light of the bifurcation theory and chaotic dynamics. 
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1.2 Literature Review 

Systems having clearances and impacts have been modelled and 
studied by several authors. It has also been understood that at 
the level of mathematical modelling, many such systems reduce to a 
finite degrees-of-freedom systems having piecewise smooth 
restoring and dissipative force characteristics. A vast amount of 
literature exists on mechanical systems with clearances and 
impacts as well as on the qualitative study of dynamics of general 
piecewise smooth systems. Nevertheless, complexities associated 
with these systems are not completely unfolded till date. Research 
is still going on and the literature is expanding. The available 
literature is broadly categorised into two parts, namely 

(i) Dynamics of mechanical systems with clearances, impacts and 
piecewise smooth non-linearity. 

(ii) Use of clearances and impacts for vibration control: impact 
dampers . 

1.2.1 Dynamics of Mechanical Systems with Clearances, 

Impacts and Piecewise Smooth Non-linearity. 

The dynamics of piecewise smooth systems has been 
studied mainly through single and multi degree (s) -of-freedom model 
having bilinear/trilinear restoring and damping force 
characteristics. Some amount of theoretical results on the free 
oscillations of piecewise linear systems can be found in the book 
by Minorsky [1] under the heading of discontinuous theory. Maezawa 
et al. [2-4] have discussed different types of non-linear 
resonances of a periodically forced, single degree-of-freedom 
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oscillator with piecewise linear restoring force characteristics. 
Watanbe [5] analysed a single degree-of-freedom piecewise linear 
system having triangular hysteresis loop. Tomlinson et al. [6] 
have studied the frequency response characteristics of structures 
with single and multiple clearance type non-linearity. Shaw et 
al.[7], Natsiavas [8-10] and several other researchers [11-15] 
have thoroughly investigated different types of non-linear motion 
including chaos in bilinear forced oscillators. A detail account 
of the complex motions of offshore structures modelled by 
piecewise linear oscillators can be found in reference [16]. It 
has been demonstrated that even a simple system like a bilinear 
oscillator can show superharmonic and subharmonic resonances and 
different types of instabilities leading to chaos. Shaw et al.[17] 
have discussed the dynamic response of a preloaded joint. Choi et 
al.[18] considered the non-linear response of an articulated 
loading platform (ALP) tower with the help of a piecewise 
non-linear, single degree-of-freedom model. The dynamics of a beam 
with motion limiting elastic constraints has been analysed both 
theoretically and experimentally by several researchers [19-20]. 

A good deal of work has been done addressing the frequency 
response characteristics of two colliding masses known as an 
impact pair. Comparin et al.[21] have studied the spectral 
characteristics of an impact pair under harmonic excitation. In 
this investigation, the impact between two masses has been 
modelled by a spring-dashpot which eventually gives rise to a 
single degree-of-freedom system having piecewise linear restoring 
force characteristics. Bapat et al. [22-23] addressed the same 
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problem where the impact is modelled by a rule (using the concept 
of coefficient of restitution) governing the instantaneous change 
in the velocity after the collision. The dynamics of an impact 
pair under random excitations has also been investigated [24-26]. 
The dynamics of an inclined impact pair [27] has shown some 
fundamental differences from that of a horizontal one. The global 
dynamics, stability and chaos of an impact pair have been studied 
by Holmes [28] with the help of a non-linear map. The same problem 
is also treated by Bapat et al.[29] and others [30]. 

When one of the masses involved in an impact pair possess 
finite stiffness and inertia while the other has infinite 
stiffness and inertia the system is known as an vibroimpacting 
oscillator. The motivation behind the study of such systems stems 
from the practical need of designing impacting tools. Senator et 
al. [31] and Fu et al. [32] investigated the noise stabilization 
problem of a vibrating hammer, Whiston [33] and Fang et al.[34] 
considered a single degree-of-f reedom impact oscillator under 
periodic excitation. Nguyen et al. [35-36] have carried out a 
detail parametric study and discussed some non-dimensional design 
parameters for the impact behaviour of an oscillator with limiting 
stops. Whiston [37] has considered the dynamic singularities 
associated with the global dynamics of an impact oscillator. 
Experimental investigation on the dynamics of an impacting 
pendulum system has been carried out by Moore et al. and others 
[38-39]. Aidanpaa et al. [40] considered a two degrees-of-freedom 
linear oscillator, impacting against a rigid stop, to investigate 
the periodic and chaotic responses. Recently a new type of 
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bifurcation has been detected in an impacting oscillator by 
Nordmark [41] and others [42]. The chaotic motions of impacting 
systems have also been addressed by Thompson et al.[43] and 
Peterka [44]. Hendricks [45] considered the dynamics of an impact 
print hammer of a dot matrix printer and demonstrated the 
possibility of the occurrence of chaos impairing high frequency 
printing. Several investigations have been carried out on the 
dynamics of impacting beams [46-48] and strings [49]. Jing et 
al.[50] reported the random dynamics of a single degree-of -freedom 
vibroimpacting oscillator where the impact has been modelled by a 
spring having Hertzian type contact stiffness. Palej et al.[51] 
have taken a different approach for analysing the impact 
oscillator where the impact force is treated as a sequence of 
Dirac distribution. Zhuravely [52] used the method of non-smooth 
transformation for analysing such a problem. Davies [53] 
investigated the stochastic vibration of an impacting beam. Bapat 
et al. [54] have recently shown the dynamic similarities among 
several vibroimpacting systems. 

A general treatment of the dynamics of linear multi 
degrees-of-freedom system with colliding components can be found 
in a paper by Natsiavas [55]. In this paper, some problems have 
also been solved with special application to vibration control. 
Padmanabhan et al. [56] resolved some spectral coupling issues 
associated with a two degrees-of-freedom system with clearance 
type non-linearity. The frequency response characteristics of 
multi degrees-of-freedom system with clearances has been studied 
by Comparin et al.[57]. The non-linear dynamics of a rotor with 
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bearing clearance is reported in reference [58]. Kaharaman and 
Singh [59-60] have studied the non-linear dynamics of geared 
systems in the presence of time varying mesh stiffness and 
clearance. The effects of multiple clearances in a geared rotor 
bearing system have also been investigated by the same authors 
[61]. A theoretical and experimental study of the dynamics of a 
four bar chain with bearing clearances is reported by Stammers et 
al.[62]. The chaotic dynamics of heat exchanger tubes in the 
presence of support clearance has been investigated by Paidoussis 
[63-64] . 

1.2.2 Use of Clearances and Impacts in Vibration Control; 

Impact Dampers 

Analytical and experimental studies have demonstrated 
that a class of non-linear auxiliary mass dampers, classically 
known as impact dampers, can be effectively used to neutralise the 
vibrations of various systems like light poles, turbine blades, 
antennae, aircraft wings, chimneys, machine tools etc. It is known 
that a suitable non-linear vibration neutraliser is, in general, 
more effective than a linear vibration neutraliser. But, in 
practice, as mentioned earlier, the introduction of a suitable 
non-linearity in the system may be difficult. In an impact damper, 
clearance and impact non-linearity is used which is comparatively 
easier to generate and control. The idea of using impact dampers 
for vibration control was conceived as early as late forties of 
this century. Since then, numerous papers have been published 
discussing its design and performance. Initially, impact dampers 



were called acceleration dampers [65]. The first systematic 
analytical basis for the design of an impact damper is given by 
Grubin [66] and which was later simplified by Warburton [67]. 
Masri et al. [68-69] developed an efficient methodology for 
studying the general motion and stability of an impact damper. 
Their methodology in fact has provided a beginning for the deeper 
analysis and estimation of the performance of impact dampers. 
Impact dampers have been successfully used in attenuating the 
vibrations of steam turbine buckets [70], tall flexible structures 
[71] and machine tools [72] etc. A good account of the design 
methods hnd constructional features of the impact dampers can be 
found in the book by Mallik [73]. The effect of gravity on the 
performance of an impact damper has been studied by Sadek et 
al.[74]. The stable performance of an impact damper heavily 
depends on the suitable choice of the design parameters such as 
the clearance, the mass ratio and the coefficient of restitution 
of the impact pair etc. Dokinish et al.[75] discussed the optimum 
design parameters for an impact damper. Pinotti et al.[76] 
provided design procedure and charts. Bapat et al.[77] carried out 
a detail parametric study and reported the effects of different 
design and system parameters on the performance of a single unit 
impact damper towards controlling the forced and and free 
vibrations of a single degree-of-freedom linear oscillator. 
Mansour [78] considered the effect of friction on the performance 
of an impact damper. Yasuda et al. [79] have provided an estimate 
of the equivalent damping effect of an impact damper during free 
vibration. 
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All the above mentioned investigations were performed 
assuming the excitation to be strictly harmonic which does not 
represent reality. In reality, the excitation is normally random 
in nature. Masri et al. [80] considered the response of an impact 
damper under stochastic excitations. Recently Semercigil et 
al.[81] have also treated the same problem. 

Despite several attractive merits viz, simple constructional 
features, easy tunability and high amplitude reduction capability 
etc. the impact dampers have got certain limitations. The optimum 
performance of an impact damper is impaired in the presence of 
broad band excitations. Several attempts have been made to 
increase the operational bandwidth of an impact damper. One 
attractive way is to use a conventional dynamic vibration 
neutraliser where the auxiliary mass is provided with motion 
limiting stops [82-83]. Such an arrangement makes use of the good 
features of both the conventional dynamic neutraliser (larger 
operational bandwidth) and the impact damper (high attenuation) . 
The use of spring supported secondary mass is another alternative 
[84]. Recently, Semercigil et al.[85] proposed a new vibration 
absorber which uses a conventional tuned absorber and an impact 
damper. Such a device is demonstrated to offer better vibration 
attenuation characteristics under broad band excitations. 

The quest for quieter technologies gave birth to bean bag 
damper [86-87] where instead of a rigid slug, packaged particles 
known as bean bag is used as the secondary loose mass. Due to the 
plastic nature of the bean bag, the noise produced by an impact 
gets reduced. Recently the concept of migration control using an 
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impact damper for a non-linear primary system is proposed by 
Chatter jee et al.[88]. Such a scheme obviates the need of 
continuous noisy operation of an impact damper. 

Another way of enhancing the performance of an impact damper 
is to use multi-particle single-unit or multi-unit systems. The 
motion and stability of an impact damper consisting of 
two-particle single unit container, have been discussed by Masri 
et al. [89]. Analytical and experimental investigations of 
multi-unit impact dampers have been reported by several 
researchers [90-92]. Cempel [93-94] used the continuous force 
approach and receptance model for analysing the performance of a 
multi-unit impact damper (MUVIN) . 

Most of the above investigations were carried out for a 
single degree-of-freedom model of the primary system which clearly 
signifies the importance of a single dominant mode of vibration. 
But such an assumption is true only for a simple system with well 
separated modes. When several modes are of significance, a multi 
degrees-of-freedom model of the primary system is necessary. 
Keeping that in view, the potential use of an impact damper for a 
multi degrees-of-freedom system has been studied by Nigam et 
al.[95]. They have questioned certain basic assumptions made by 
the earlier authors. Recently, Chelmers et al.[96] used two-unit 
impact dampers, placed at different locations, for controlling the 
vibrations up to the second mode of a cantilever beam. 

The use of motion limiting stops for a centrifugal pendulum 
vibration absorber has also been reported [97]. 
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Despite several modifications, a passive impact damper is 
incapable of performing well in a varied dynamic environment. This 
fact led Masri et al. [98-99] to develop an active parameter 
control technique for controlling the vibrations of an arbitrary 
non-linear system under a general kind of loading. One of the 
strong features of the technique is that no mathematical model of 
the system is at all necessary to design the control strategy. 

It is a well known fact that even a low dimensional model of 
a non-linear system can very well explain the existence of 
different kinds of complex instabilities leading to chaos. On one 
hand the study of chaos gives some explanations of the presence of 
noise in an otherwise simple system and on the other hand, it 
gives the indications of potential dangerous components lying 
hidden in the beneficiary applications of non-linearity. Shaw et 
al.[100] have studied the nature of chaotic motions of a two 
degrees-of-freedom impact oscillator. Recently Sung et al.[101] 
have reported the existence of chaos in the model of an impact 
damper . 


1.3 Objectives and Organization of the Present Work 
Very often vibratory systems are modelled as linear systems. 
The principles and hardware designs of the control methods for 
such linear vibrations are well documented. A design based on the 
linear model suffices if the oscillatory movement is small. 
However, the non-linearity becomes pronounced when the vibration 
is large and needs to be controlled. Furthermore, there exist real 
life systems which, due to the presence of clearances and 
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self-excitation mechanisms, are basically non-linear. The 
inclusion of a secondary system in the form of a loose mass which 
undergoes repeated collisions with the primary system is a common 
method of vibration control. Such a secondary system is known as 
an impact damper. The major objective of the present thesis is to 
study the performance and design of an impact damper attached to 
non-linear primary systems. Three different models of the primary 
system, viz., 

(i) harmonically excited Duffing^s oscillator, 

(ii) autonomous van der Pol oscillator 

and (iii) harmonically excited van der Pol oscillator 
have been investigated. Both elastic and inelastic impacts have 
been discussed. Besides the design of an impact damper for these 
systems, the bifurcations and chaotic dynamics of the last two 
systems have been given special attention. Accordingly, the 
presentation has been organized as detailed below. 

The systems under consideration can be classified as 
piecewise smooth, non-linear systems. In chapter 2, the equivalent 
linearization procedure has been extended for such systems and is 
termed as the equivalent piecewise linearization (EPL) method. The 
efficacy of the method, as compared simple harmonic balance method 
(HBM) , for piecewise smooth, non-linear systems has been 
demonstrated by considering a Buffing's oscillator with motion 
limiting stops (Vibroimpacting oscillator) . The stability analysis 
of the solutions is also presented. All the analytical results are 
verified by numerically integrating the equations of motion. 


In chapter 3 , the design and performance of an impact damper 
for controlling the high-amplitude vibration of a harmonically 
excited, duffing type primary system are discussed. Again, the 
results are obtained by HBM, EPL and numerical integration. A 
detailed parametric study is carried out to demonstrate the 
effects of various parameters. The procedure for arriving at an 
optimal design is presented. The question of controlling the 
subharmonic response has also been addressed to. Finally, the 
possibility of using the impact damper discontinuously in an 
active manner has been demonstrated. In this mode, the impact 
damper is shown to act as a migration controller which changes the 
initial conditions so as to bring the original system within the 
basin of attraction of a lower amplitude solution. Thus the 
feature of multiple stable solutions of a non-linear primary 
system is put to use. 

In chapter 4, the performances of an impact damper towards 
controlling the autonomous (limit cycle amplitude) and 
harmonically entrained forced oscillations of self-excited 
oscillators, are discussed separately. In the latter case, after a 
detailed parametric study, an empirical relation for the optimum 
design is established. A sensitivity analysis of the design with 
reference to the self-excitation and forcing parameters is 
included. 

In chapter 5, a study of bifurcations and chaos in 
self-excited, autonomous oscillators with an impact damper is 
presented. Results have been obtained by considering various 
models of the self-excitation mechanism. The major emphasis is on 
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the detection of the types of bifurcations which are independent 
of the exact nature of the model. In this chapter, a numerical 
technique for analysing the eigen properties of the numerically 
integrated response has been proposed. 

In chapter 6, different routes to chaos for a forced, van der 
Pol oscillator with an impact damper have been investigated. Some 
uncommon intermittency routes to chaos have been detected. Finally 
in chapter 7, the overall conclusions and some suggestions for 
future works are included. 



CHAPTER-2 


PERIODIC RESPONSE OF PIECEWISE SMOOTH NON-LINEAR 
OSCILLATORS UNDER HARMONIC EXCITATION 

2.1 Introduction 

A harmonically excited non-linear oscillator has been 
studied extensively in various branches of mathematics, physics 
and engineering. More commonly, the non-linearity in such 
oscillators is represented by the restoring and/or damping forces 
which are smooth non-linear functions of the state variables. 
Another source of non-linearity is the presence of discontinuity 
in the stiffness or damping coefficient or in the state variables. 
There exist a number of real life oscillators whose mathematical 
models include both the sources of non-linearity mentioned above. 
Typical examples in this class include mechanical devices with 
clearance and backlash, non-linear oscillators with 
motion-limiting stops or impact dampers, beams with cracks etc. 
Such an oscillator, for which the restoring (and/or damping) 
force is represented by a number of smooth non-linear functions, 
each valid for a range of the state variables, is called a 
piecewise smooth, non-linear oscillator. The harmonic balance 
method (HBM) is most commonly used for obtaining the periodic 
solutions of a harmonically excited, non-linear oscillator. It has 
also been used for piecewise smooth, non-linear oscillators [88]. 
The HBM, though works quite well for a weakly non-linear system, 
fails to give satisfactory results if the effective non-linearity 


is large due to the sharpness of the change between the adjacent 
smooth functions. Moreover, the inclusion of higher harmonics in 
closed form becomes impossible due to the presence of Heaviside 
step functions in the mathematical model of such an oscillator. 

The difficulties encountered in the HBM while tackling a 
piecewise smooth, non-linear oscillator are taken care of by the 
incremental harmonic balance method (IHBM) [102]. This method, 
based on the local linearization in the parameter space, is also 
capable of including higher harmonics in the solution. Kim et al. 
[11] have developed an useful method by combining the fast Fourier 
transform and harmonic balance techniques (FFT-HB) . This method 
has been used to study the dynamics of both piecewise linear and 
piecewise smooth, non-linear systems [18]. 

The accuracy of the results obtained by the IHBM and FFT-HB 
technique depends, like all harmonic balance techniques, on the 
number of harmonics included. Sometimes, ten to twenty harmonics 
are to be retained for convergence [18]. Basically both the 
methods are iterative and hence computationally involved. No 
uniform convergence criterion can be ‘used in the iterative scheme. 
These problems become more prominent if the response contains 
sharp spikes (e.g., due to a sharp change in the restoring force 
as in the impacting situation ) when a large number of harmonics 
need to be retained. The stability analysis associated with the 
above methods uses variational technique for non-dif ferentiable 
systems (which is mathematically heuristic) and gives rise to 
Hill's equation. Both the IHBM and FFT-HB technique, while having 
general applicability, may turn out to be computationally 
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inefficient for some simple non-linear systems undergoing sharp 
changes in the restoring and/or damping force. Sometimes recasting 
of a problem to suit the algorithm also becomes complicated. 

In the present chapter, we propose to first transform a 
piecewise smooth, non-linear system to an equivalent piecewise 
linear (EPL) system and thereafter, use an exact method proposed 
by Masri el al. [69]. In this method, the form of the solution for 
each range having different coefficients of equivalent linearity 
is known. The solution contains a number of unknowns which are 
subsequently determined by imposing the conditions of periodicity, 
symmetry (if any) and consistency to be maintained at every 
junction of adjacent linear ranges. Thus, in the present method, 
the form of the solution is not based entirely on the harmonic 
components and the sharp change between the adjacent non-linear 
ranges are treated as such. Consequently, the sharp spikes in the 
time series of the responses are predicted much more easily as 
compared to that by the IHBM or FFT-HB technique. A dual 
representation of the solution is used to greatly enhance the 
algebraic simplicity and numerical efficiency. The stability 
analysis is carried out using the method of error propagation 
[68] wherein a linearized Poincare' map relating the perturbation 
vectors at the end and beginning of a period is constructed. The 
eigenvalues of the Poincare map determine the stability of the 
obtained solution. 

It must, however, be mentioned that an accurate knowledge 
of the number of junctions to be encountered during one period (of 
the periodic response) is a pre-requisite for using the proposed 
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method. The method obviously becomes unwieldy if the number of 
such junctions becomes too many. Fortunately, for a number of 
models of practical interest, this number is often small for a 
wide range of parameter values. Numerical results are included to 
demonstrate the efficacy of the proposed method for the following 
two situations; 

i) a harmonically driven Duffing' s ossillator having an 
one-sided elastic constraint and 

ii) a harmonically driven buffing's oscillator having an one 
sided rigid constraint. 

2.2 Harmonically Driven Buffing's Oscillator with One-sided 
Elastic Constraint. 

To demonstrate the efficacy of the proposed method (EPL) , 
let us consider a harmonically excited Duffing' s oscillator with 
one-sided elastic constraint (Figure 2.1). The equation of motion 
of the system in non-dimensional form is given by 

X+2hX+X+EX^+ r^{X-l) U(X-l) = F^cos nr (2.1) 

where X = ^/^q/ K^/m, h = C / ( 2 m , 

= = K2>'o ' “o> ' = K3 / (m 0,2) , X = o,^t, 

n = 01/0,^, = P /(m 0,2 Xjj) 


and the dot denotes differentiation with respect to the 




oscillator with an one-sided elastic constraint. 
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non-dimensional time t. The Heaviside's step function U(.) is 
defined as 


U(Z) =1 Z > 0 

= 0 Z < 0 

It is desired to suitably transform equation (2.1) into piecewise 
linear form given by 

• « « 

X-2pX+X+eKgX+r^ (X-1) U (X-1) = cos Qt (2.2) 

With p = - h and as an unknown constant. Equation (2.2) can be 

solved using the method of condition matching [69]. However, the 

solution so constructed carries the unknown K , which can 

e 

subsequently be calculated as elaborated in the next section. 

2.2.1 Periodic Solution 

We are interested in constructing a periodic solution of 

equation (2.2) having a period r , which is an integral multiple 

P 

of the forcing period (x^ = 2uJQ.) . To start with, we consider the 
type of solution depicted in Figure 2.2. Numerical simulation of 
equation (2.1) suggests that this type of motion of the oscillator 
does occur for a wide range of parameter values. In Figure 2.2, we 
have arbitrarily set x = 0 when the contact with the constraint is 
first established with positive velocity and the contact is 
maintained for a time interval x^. This arbitrary setting of the 
time origin requires introduction of an unknown phase (p in the 
forcing function. Thus, a solution of the form shown in Figure 2.2 
can be written as 
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X (T) = X^(T) U (-T) + X2(t) U (T) (2.3) 

where T = t - and X^ and X^ are given by (see equation (2,2)) 


’ ■ * 2 

X^ - 2p X^ + ( 0 ^ X^ = cos (Qr + <p) + (2.4) 

with X (t) > 1 , -¥- 0 < T < T 
X c 

and 

X^ - 2p X 2 + ^2 X 2 = cos (Qt + (p) (2.5) 

with X (T)<1, ■V0<T<T<T 

o p 

with and defined as 


cj^ = 1 + eK^ + r. 

1 el 

wj = 1 + eK 

2 e 

The solutions of equations (2.4) and (2.5) can be written as 


^ (t) = e^^ (d^ cos + d 2 sin(‘r}^T)) + M^cos(nT + (p + 0^) 




( 2 . 6 ) 


and 

X 2 (T) = e (d^cos t} 2 (t-t^) + d^sin 

+ M2COs(nT+ <^+02) (2.7) 



where 
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M. = FJ / (w? - + (2pn)^ 

and tan 0. = ^ ^ with i = 1,2. 

^ w. - 
1 

Satisfying the continuity condition (at x = r ) and using the 
periodicity with time-period one can write 

Xi(0) = (x^) = X^ (x^) = X^ (Xp) = 1 

Xi(0) = X 2 (Xp) = X^ (say) 

and X^ (x^) = X^ (x^) = X^ (say). (2.8) 

It should be noted that equations (2.6) and (2.7) contain seven 
unknowns, vi 2 ., d^, d^, d^ , <j>, x^ and K^. Using conditions 

(2.8) we can eliminate d^^, d^, d^ , d^ and 0 to obtain a relation 
of the form 

ff (K^, X^) = 0 . (2.9) 

One may think of several ways of computing from (2.9) with the 
objective of minimizing the integral 

I = J X(x) “ 

o 


2p X(x) + X(x) + c X (x) + 

X 2 

rj^(X(x) - 1) U(X(x) -1) - Fq cos(nx) dx 
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where X(t:) is defined by equation (2.3). However, direct 
substitution of equations (2.6) and (2.7) in the integral I makes 
the subsequent algebra too formidable. An iterative way would be 
to start with an initial guess value of K^. But such a numerical 
procedure is found to be computationally very inefficient simply 
because no good initial guess of can be made in a systematic 
way and the solution of equation (2.9) along with the evaluation 
of the integral I makes the algorithm very slow. 

We, therefore, suggest the following method. First the 
solution X(t:) is expanded in its Fourier basis and written as 

X*(t) = cos (Qt + i//^) + (2.10) 

and is determined as a function of the Fourier coefficients 

^e = ^e (^o' ■^1' ) (2.11) 


to render 



K X* )2 dT 
e ' 


0 . 


( 2 . 12 ) 


It should be noted that X (x) and X(x) are two different 
representations of the solution of the same differential equation. 
So these should be brought closer by minimizing the least square 
of their difference, i.e.. 
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Tp 

J (X*(t) - X(r))^ dr is minimum (2.13) 

o 

According to the theory of Fourier series, this minimization takes 

place if A^, A^..etc. correspond to the Fourier coefficients of 

X(t) given by equation (2.3). This condition generates appropriate 

number of equations to solve the problem as explained below. 

Let us consider 

X* (r) = A^ + A^ cos (fix + (2.14) 

when equation (2.12) yields 


where 


and 



2f A + f, A 
a o b 

2 2 

2A + a; 

O 1 


3 3 2 

% + I ^o ^1 


2 3 

3 A A, + T A 
o 1 4 


1 


3 

1 • 


From condition (2.13) we get 



2 


(2 . 15) 


(2.16) 


and 


. dx 


(2.17) 
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where the constants appearing in X(t) (see equations (2.6) and 

(2.7)) are also expressed in terms of K and t . Equation (2.9) 

6 o 

can be rewritten, by using equation (2.15), in the form 

f2 (A^, A^, = 0 (2.18) 

Now equations (2 . 16) - (2 . 18) can be solved to obtain A^, and 
when all the quantities appearing in equations (2.6) and (2.7) can 
easily be obtained. 

It should be emphasized that in this method one gets dual 
approximations of the solution; namely X (t) and X(t) . So far as 
the nature of the solution is concerned, X(r)^ gives better 
approximation as compared to X (t) . However, it will be seen later 
that sometimes X (r) can also accurately predict the maximum 
displacement of the oscillator. We call X(r) as the primary 

k 

solution and X (r) as the secondary solution. 


2.2.2 Stability Analysis 

Both stable and unstable solutions are obtained by the 
proposed EPL method. The method of error propagation explained by 
Masri et al. [68] has been used to ascertain the stability of the 
obtained solution. Towards this end, we consider the obtained 
solution X(t) represented by the solid line in Figure 2.3 and its 
perturbed form X(t) represented by the dashed line in the same 
figure. Referring to Figure 2.3 we can write 
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X^(0) = (T^ + AT^) = 1 

• • • 

X, (O) = X + AX 
1 ' ' o o 

• • * 

X, (T + AT ) = X + AX (2.19) 

1 ' c o' c c ' ^ 

• 9 

(see condition (2.8) for the definitions of X and X ) 

' ' o c 

where AT = Ar - Ax (2.20) 

o c o 

Using the above conditions one obtains the linearized relations 
between the perturbations in the form 


f ] 


f ^ 

Ar 


Ax 

c 

' = [P.] ^ 

0 

Me 


Me 

c 

V J 


0 

\. / 


where the elements of the 2x2 matrix functions of the 

system parameters and the constants appearing in the solution 

given by equation (2.6). Considering the solution for the interval 

beyond t (and ^ x ) given by equation (2.7) and using the 

c p 

following conditions; 

X2(0) = X 2 (Xp + AT^) = 1 , 

^2(0) = 

and X 2 (Tp + ATj^) = X^ + AX^ (2.22) 

where AT, = Ax - Ax - x and X- is the velocity of the oscillator 

1 p c c f 

at the end of the period, one obtains 
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fix 1 


r \ 

Ax 

> 


c 

> 

Ak 

c 


Combining equations (2.22) and (2.23) one gets 


r 1 

Ar 


' \ 

At 

p 

' = [P] ^ 

o 

V 



Ak 

f 

v .y 


Q 

/ 


where [P] = [P^] [P^] . A solution is stable if all the eigenvalues 
of [P] lie inside the unit circle around the origin of the complex 
plane. 

2.2.3 Harmonic Balance Techniques 

It is known that the IHBM, FFT-HB and HBM are essentially 
same in the first approximation. As already stated, while the 
first two methods can be extended to higher order approximations, 
the HBM can include only the first term in closed form. To see the 
efficiency of the proposed EPL method (in the first term 
approximation only ) the HBM has been chosen for comparison. Let 
the first term solution in the HBM be 


X (t) = B^ + Bj^ cos (nr + 0) (2.25) 

Substituting equation (2.25) in equation (2.1) the 
coefficients of equal circular functions are matched. While doing 
so, the Fourier expansions of the Heaviside's step functions are 
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used. Finally one gets the following non-linear algebraic 
equations 


+ eB^ + ^ B^ bJ + r, (D =0 (2.26) 

0020ll0 

B^ (1 - n^) + ^ + 3 e B^ B^ + r^ 0^ = cos 8 (2.27) 

and 

Where, 
then 


-2 h B.n = F sin 8 
1 o 


(2.28) 


if the collision occurs, i.e., if 


®o - " 

B. 


1 , 



and (D^ 



otherwise 0 = Q, = 0. 

o 1 

The solution X(T) is thus obtained by solving B^ and B^ from 
equations ( 2 . 26 )-( 2 . 2 a). 

2.2.4 Numerical Results 

The periodic solutions with period 2Tr/f2 have been obtained 
using both the EPL method and simple HBM. These results are 
compared with the solutions obtained via direct numerical 
integration using fourth-order Runge-Kutta-Merson algorithm with 
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adaptive step size control. A typical frequency response 
characteristic of the system is shown in Figure 2.4. The branches 
AB and DE in Figure 2.4 correspond to the non-contacting solution 
and hence can be produced using several techniques widely 
discussed in the literature. Even a single term harmonic balance 
can also produce these branches accurately. The branch DC refers 
to the unstable solution. So we are interested in obtaining only 
the branch BC where the elastic constraint is in operation. 

The response of the system being asymmetric, the maximum 
displacements of the mass on either side of the equilibrium 
position (i.e. both positive and negative) are plotted separately 
in Figures 2.5(a) - 2.5'(d) for two different values of r^. 

The maximum displacement on the free side is represented by 
5^ as shown in Figure 2.2 whereas that on the constrained side is 

"k 

S^. It may be pointed out that the secondary solution X (t) gives 
another measure of the maximum displacement of the mass on the 
free side by S* = I X* I max' 

From Figures 2.5(a) and 2.5(b) it may be seen that for a 
low value of the stiffness of the restraining spring / both 
the EPL and HBM produce a good estimate of 5^ which is also seen 
to be higher than 5 . The EPL method, when compared to the HBM, 
gives better results so far as the displacement on the constraint 
side is concerned. This will be of significance especially if the 
deformation of the restraining spring is an important parameter 
for its design. From Figures 2.5(a) and 2.5(c) we note that as the 

•k 

displacement on the free side increases, 5^ becomes more accurate 
than 5^. From Figures 2.5(c) and 2.5(d), we can conclude that with 




Figure 2.5 Frequency response within the contacting branch near 

fundamental resonance, e = 0.05, F = 1, h = 0.05 r = 

o 1 

5 in a & b and r*- 25 in c & d . EPL, - - - HBM,- 

~ - 5 o numerical integration. 
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a stiff restraining spring, the single term EPL solution is much 
better than that produced by the HBM. The accuracy of the HBM was 
found to deteriorate drastically with increasing values of r^. 

2.3 Vibroimpacting Buffing’s Oscillator 

In the previous section we have discussed the efficacy of 
the EPL method in obtaining the approximate solution of a 
Buffing's oscillator with one sided elastic constraint. In the 
present section, the constraint is assumed to be rigid. A rigid 
constraint can be modelled in various ways. One way would be to 
assign a very high value (ideally infinity) of r^. But a very high 
value of r^ leads to unavoidable numerical difficulty in solving 
the algebraic equations (2 . 16) - (2 . 18) . This difficulty is overcome 
if the collision of the oscillating mass with the rigid constraint 
is modelled by assigning a rule governing an instantaneous change 
in the velocity (using the concept of coefficient of restitution) 
at the constraint boundary. It should be noted here that such 
modelling is not amenable to the HBM. However, the EPL method can 
be used very efficiently as explained below. 

2.3.1 Periodic Solution 

For the sake of algebraic simplicity, first a single impact 
per cycle motion having a period t is constructed. The nature of 
the solution sought is sketched in Figure 2,6 when the 
displacement of the mass can be written as 
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X(t) = (d^ cos(T} 2 r) + sin {-q^x)) + cos (CIt + (p + 6^) 

(2 . 29) 

Further, we have the following periodicity and impact conditions: 

X (0) = X (Tp) = 1 

let X (O"^) s X and X {x~) = X. (2.30) 

' a ' p D 

with ia = - * ib 

where x is the coefficient of restitution. 

Using conditions (2.30) we eliminate two unknowns, viz.,d 2 and 
to obtain a relation of the form 

K^) = 0 (2.31) 

For calculating K corresponding to a solution of period t = 2nlQ, 

e p 

we consider the secondary solution 

X*(T) = cos (nr + 5^) + A 2 cos (2nT + 3^) (2.32) 

when equation (2.12) yields 

“a *0 + *d ''2 + <*1 <^2 ' 

.. 2 , , 2 ,, 2 
2 A- + 

0 12 
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Where f, = J A? A + J A^ A + ~ A^ 

a o2lo 22o 4 2 


„ 3A^ A 

fj^ = 3 a; A^ + I A^ + — 


f_ = 


3 A^ A^ A^, 




f = 


3 ,2 ^ 

^ A. A^ 
2 10 


and equation (2.13) yields 



and 


tan 6 

n 


Im (X^) 
Re (X^) 


for n = 1,2 


where 



iflnr 
0 «> 


dr 


and Im (.) and Re(.) represent, respectively, 
real parts of X^. 


A^ cos (§2 - 25^) , 


(2.33) 


(2.34) 


(2.35) 


the imaginary and 
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A similar procedure can be followed for obtaining a 
solution with higher (than that of the excitation) time-periods. 
For example, to compute a solution with period 2 t^ the secondary 
solution is written in the form 

X*(T) = + A^ cos {Qx + 5^) + A^^^cos (^x + 5^^^) 
when finally we obtain 


K 


= { 


f A + 
a o 


^b 




(A, 


+ 


^l/2^h^ 


cos (5. 




D 


where t = A^ + ^ A A^ + ^ A^ . A + |- A^,„ A. cos (5 - 25. 

a o 2 o 1 2 1/2 o 4 1/2 1 2 1/2 




and 


- _ 2 7\2 T, 

^c 2 ^1/2 o' 


4 ^1/2 ^^o ^1/2 2 ^1/2 ^1 ' 


f^ = 


3A„ A, 
o 1 


D = 2 aJ + A? + A?,- 
O 1 1/2 


The stability analysis of the solution is carried out using 
the method of error propagation outlined in section 2.2.2. 
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2.3.2 Numerical Results 

In this section, the numerical results obtained from the EPL 
method are compared with those obtained from direct numerical 
integration. Only the stable impacting solutions are plotted. The 
maximum displacement of the mass on the free side, (5^ and 5*}, 
obtained with t = 2TT/n is plotted against frequency in Figure 
2.7. It is seen that the EPL approximation with 0 is better 

than that with = 0. In fact, the inclusion of more number of 
harmonics in X (t) is seen to produce higher accuracy. However, 5^ 
obtained with A 2 = 0 gives a reasonably good estimate of the 
maximum absolute displacement of the oscillation, especially if 
the value 

* * 

of is high. The impact velocity (Xj^) is plotted against the 

forcing frequency in Figure 2.8. It is seen that even the first 
approximation (i.e. with A 2 = 0) gives a very accurate estimate of 
the impact velocity. 

The maximum absolute displacement 5^ for a solution with 
period An/Q is plotted in Figure 2.9. Again it is confirmed that 
the EPL method can very well account for the subharmonic 
contributions . 


2.4 Conclusions 

Simple single term harmonic balance method fails to predict 
the periodic response of a piecewise smooth, non-linear 
oscillator, especially if the degree of non-linearity changes 
sharply at the junction between the adjacent (non-linear) regions. 
The proposed equivalent piecewise linearization technique works 



Figure 2.7 Maximum displacement of the vibroimpacting Duffing’s 

oscillator with x = 2n/Q. x = I values of the 

P 

parameters h, F'^and c are specified in Figure 2.5. — 

- - 5 ^ (A^ = 0), ^5^ (A^ * 0), - . - . - sj (A^ = 

0), o numerical integration. 




4 



Figure 2.8 Variation of the impact velocity with forcing 
frequency. Values of h, F , c and x are specified in 

Figures 2.5 and 2.7. EPL (A^ = 0), o numerical 

integration. 




Figure 2.9 Maximum displacement of the vibroimpacting Duffing’s 

oscillator with t = 4t[/Q.. Values of h, and c 

p o 

are specified in Figures 2.5 and 2.7. o 

numerical integration. 
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much better for such problems. The main reason for this can be 
explained as follows. The success of the HBM, being basically a 
frequency domain technique, depends heavily on the inclusion of 
all the prominent frequencies in the assumed form of the solution. 
However, for the present problem this is not possible while using 
the simple HBM. We may add that there may be situations when the 
simple HBM may provide a good estimate of the maximum response 
while completely failing to predict the time response. Such an 
example will be discussed in the next chapter. The proposed EPL 
method, on the other hand, simultaneously utilizes both the time 
domain and frequency domain information. This method not only 
gives better results than that obtained from the simple harmonic 
balance method but is also capable of including super and 
subharmonics in the solution. A dual representation of the 
solution greatly simplifies the algebra involved in the proposed 
EPL method. The method is not iterative hence computationally less 
intensive for simpler systems. The EPL method becomes 
algebraically unwieldy only if the number of smooth non-linear 
regimes increases. The stability analysis of the obtained solution 
can be carried out by the method of error propagation. 



CHAPTER 3 


IMPACT DAMPERS FOR NON-LINEAR VIBRATING SYSTEMS 


3.1 Introduction 

In this chapter, the performance and design of an impact 
damper for controlling the high-amplitude vibration of a 
non-linear oscillator is discussed. Specifically, the primary 
system is taken as a harmonically excited, viscously damped, 
Duffing's oscillator. Two different approaches are considered. 
First, the secondary mass is used in a passive manner to generate 
a stable low-amplitude solution for the primary system. Towards 
this end, an initial design, obtained by the single term harmonic 
balance method, is improved upon using the EPL technique and 
numerically integrated results. In the second approach, the 
multistability feature of the non-linear oscillator is utilized. A 
guided "on-off" control strategy is followed to use the secondary 
mass for changing the initial state of the primary oscillator. 
This change of state brings the primary oscillator within the 
basin of attraction of the non-resonance (low amplitude) branch 
from the resonance branch. This procedure has been referred to as 
the guided branch transition. Numerical results are presented to 
verify the solutions, their stability and thereby the feasibility 
of the proposed approaches. 
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3.2 Theoretical Analysis 


3.2.1 Harmonic Balance Method 

The model under consideration is shown in Figure 3.1. 
For a linear primary system, it has been shown in previous works 
[69,77] that an impact damper undergoing elastic collisions 
performs better than the one having inelastic collisions. 
Accordingly, here also first we consider elastic impacts between 
the primary mass m^ and the secondary mass m 2 - The elastic impacts 
are modelled by two linear contact springs, each of stiffness K 2 - 
Later on, we shall include the impact damping to simulate 
inelastic collisions having a coefficient of restitution < 1. 
The value of x is known to control the frequency bandwidth of the 
effectiveness of the impact damper [69]. The friction force 
between m^ and m 2 is neglected. 

The equations of motion for m^ and m 2 are 


and 


m^x^^+ x^+ K^xJ + K 2 [(y - |) U (y 

+ (y + §) u(-y “ §) 3 ^ 

>"2y"+ K2[(y - f) U (Y - f ) 


2 ' 


(3.1) 


+ (y + f ) U(- y - f )] = m 2 X^^ (3.2) 

where the prime denotes differentiation with respect to time t, y 
= x^ - X 2 . All other symbols are explained in Figure 3.1. In 
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equations (3.1) and (3.2), U(.) is the Heaviside step function 

defined as 

U(?) ==0 if ^ < 0 and 
=1 if ? > 0. 

Substituting the following non-dimensional parameters 

^2 = ^1^0- V ' “o - , '"l = ''2/ 'Vo)' 

= m 2 /m^, d^= d/2x^, h = C^/ (m^w^) , Q = w/cj^, w = ^ 

1 / 3 , 

and = (F/K^) ' in equations (3.1) and (3.2) we get 

+ h + o)^ + r^((>(Z^) = cos Qx (3.3) 

and Z 2 + "" ^1 

where -PiZ^) = (Z 2 - d^) U (Z 2 - d^) + (Z 2 + d^) U(-Z 2 - d^) 

and the dot denotes differentiation with respect to r. 

Let us seek solutions of equations (3.3) and (3.4) in the 

form 

Z^ = a cos (Qt + 0 ^) 

and Z2 = b cos {Qt + 9^) . (3.5) 

Substituting the above solutions in equations (3.3) and (3.4) we 
expand ^ Fourier series retaining only the first 

harmonic. Thereafter, equating, respectively, the coefficients of 
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sin (Qt+ 0^) , sin (Qx+e^) , cos(QT+e^) and cos(Q'c+02) from both 
sides we get the following set of non-linear algebraic equations; 

aw^ - an^ + 3 a^/4 + r^ B(b) cos ^ ® 1 ' {3.6a) 


-h a n - r^ B(b) sin (02"®i) (3.6b) 

- b + (r^/r^^) B(b) = - a cos (9^-8^) (3.6c) 


sin (e^-e^) = 0, 


where B(b) = b 


1 


2 

TT 



,do. 2 
^ b^ 


+ 


sin 


-1 



(3 . 6d) 


= 0 



(3 . 6e) 


Equation (3.6d) implies e^- = nn, n = 0,1,2,.,.. 

If - 0^ = 0/ 211, 4n...,then equations (3.6a), (3,6b) and (3.6c) 

yield 


2 

ao)^ - 


+ 3 


a /4 + r.B(b) = cos 0. 


(3.7a) 


- h a n = sin 0^ (3.7b) 

- b + (r^/r^) B(b) = -a . 


and 


(3.7c) 
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Similarly, with 9^ - 9^ - n, Sir, 5n..., we get 

ao)^- a + 3 a^/4 - r^B(b) = - cos 9^, {3.8a) 

h a n = sin 9^ (3.8b) 

and - b + (r,/r ) B(b) = a (3.8c) 

1 m 

•ff 

which, with a= -a , can be rewritten as 

a*cj^ - a* + 3 a*^/4 + r^ B(b) = cos 9^ , (3.9a) 

- h a* n = sin 9^ (3.9b) 

and - b + (r^^/r^) B(b) = -a* (3.9c) 

Thus, we see that the solution set (a,b,9^) obtained from 
equations (3.7) is same as the solution set (a ,b, 9^) obtained 
from equations (3.9). Therefore, we use only equations (3.7) to 
get all the solutions. 

In all real life situations, the value of r^ is very 
high implying 


b = d + e , where 0 < e << 1 
o ' 

Using equation (3.10) in equation (3.6e), we obtain 


( 3 . 10 ) 
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B(b) = “ ^ ^ 1^1 (3.11) 


where e is obtained from equation (3.7c) as 


(a- d ) n' 


(3.12) 


'o 


Using relations (3.11) and (3.12) in equations (3.7a) and (3.7b) 
we finally obtain 


aw^- a n^+ 3a^/4 - r^ (a-d^) 


cos 9 . 


(3 . 13a) 


h a fJ = sin 0^ . 


(3 . 13b) 


Substituting A = tan (9^/ 2) in equations (3.13a) and (3.13b) and 
eliminating the variable a, the following equation in A is 


obtained: 


A® (u+1) + a A^ + (31/ + 1) a"^ + O + 2a) A^+(2v ~ 1) A^ 


+ a A + (f-l) = 0 


(3.14) 


where 


2 , 2 
h n 


- - Tm q" ) 


h^ 


V = r n d 
m o 


libra 

1. I. T.. KANPUB 
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After solving the polynomial equation (3.14) for A, we get the 
values of the amplitudes (a) from equation (3.13b). Before 
proceeding on to the stability analysis of these solutions 
(presented in the next section) , it may be worthwhile to note that 
equations (3.13a) and (3.13b) represent the single term harmonic 
balance solution of the following differential equation: 

* * ' 2 3 ‘ ‘ 2 

x + hx + ojnX + x = cos Qz - r x-r n d cos (nz+d. ) . 

c mmol 

(3.15) 

Equation (3.15) implies that the addition of an impact damper to 
the primary oscillator is equivalent to the addition of 
acceleration feedback and force-modification mechanisms. This 
confirms the use of the term acceleration damper [66,67] as a 
synonym for an impact damper. 

3.2.2 Stability Analysis 

To analyse the stability of the solutions obtained in 
the previous section, we consider 

= a(T) cos (fir + e^(T)) (3.16a) 

and = b(z) cos (f2r + e^(z)) (3.16b) 

with a(T), b(T), and ^ 2 ^^^ slowly varying functions of 

time. Substituting equations (3.16a) and (3.16b) in equations 
(3.3) and (3.4) and carrying out the harmonic balance after 

** »t •• •» • t 

neglecting small terms such as a(z) , b(z) , S^^Cr), 

etc, we get the following four autonomous differential equations; 
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1 / 

^ ' 2Q \ 


h a n + r, B(b) sin (0, 




+ sin0. 


(3 . 17a) 


9=- 


2 an 


{ 


cos 0 - r B(b) cos 




2 3 

4 


a^+a 


n' 


1 


(3 . 17b) 


b = - 2 ^ I sin ©2 + h a n cos 

+ (aul +1 a^) sin (©3_-02)} (3.17c) 

and 

®2 = - ^ - ’^1 + h a n sin (e^-e^) 

- (aw^+ I a^ )cos ^ 2 } * (3.17d) 

While deriving equations (3 . 17a) -(3 . 17d) , the Fourier series of 
the function (p(Z^), with = b(T) cos (fir + 02 ( 1 :)) had to be 
carried out. However, the Fourier coefficients of <P(Z^) can not be 
obtained in a closed form. Hence, the following approximation is 
made. Since r^ is very high, the contact time during an impact is 
very small. Consequently, <Pi'^2^ zero for most of the time 

during the period In/Q and assumes nonzero value only for a very 
small interval of time during which the slowly varying functions 
b(T) and 0 „(t) can effectively be considered as constants. This 
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approximation makes the closed form evaluation of the Fourier 
coefficients possible. The stability of the solutions can thus be 
judged from the signs of the real part of the eigenvalues of the 
Jacobian of the right-hand-side of equations (3.17) evaluated at 
(a, b, 0^, details of the Jacobian are given in the 

Appendix A. If the real part of all the eigenvalues are negative, 
the solution is stable and otherwise unstable. 

3.3 Design of an Elastic Impact Damper 
By design of the impact damper, we mean the determination of 
the values of the parameters r^^ and d^ for which the response of 
the primary system, i.e. the value of the amplitude (a), is 
rendered small. Before going into the details of this design, the 
following points should be mentioned: 

(i) The variation of the amplitude of the primary system with the 
forcing frequency, in the absence of the impact absorber, has the 
general trend shown in Figure 3.2 (for not too high values of h 
and . The frequency corresponding to the maximum value of the 
amplitude is defined as the peak frequency Q . 

hr 

(ii) The relative amount of non-linear contribution in the 

restoring force of the primary system can be controlled by the 
parameter . A higher value of implies a lower relative 

non-linearity and vice versa. 

(iii) To numerically integrate the equations of motion, the fourth 
order Runge-Kutta-Merson algorithm of NAG with adaptive step-size 
control was used. We have randomly selected a large number of 
initial conditions within a ball of radius 2 around the origin of 
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Figure 3.2 Typical f rSquency response curve of a sinusoidally 

driven Buffing’s oscillator. stable branch, - - 

- — unstable branch, A - one-third subharmonic, = 
0, h = 0.215 
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the phase-space. All these initial conditions resulted in motions 
where collisions occurred. There may exist certain special 
combinations of initial conditions (depending on the values of 
other parameters) where no collision occurs. In this work, no 
further attention has been given to such non-colliding motions. 

(iv) The steady state of the motion has been ascertained by 
considering the Poincare section (stroboscopic map with 2TT/n as 
the strobe period) of the phase trajectories. When the iterates of 
the Poincare map reach a fixed point, the corresponding motion 
obviously has a period 2Tt/Q. Similarly, the existence of n fixed 
points in the iterates implies a motion having a period 2Tin/Q. 
This fact is used to locate the subharmonic solutions. 

From the approximate solution given by equation (3.14), it is 
easy to see that if r^ and d^ are chosen to satisfy the relation 

r d = 1 , (3.18) 

then one of the values for the amplitude (a) turns out to be zero. 
Of course under such a situation i.e., when equation (3.18) holds 
good, the solution of equations (3.7) does not give the amplitude 
as exactly zero, but does reveal a very small amplitude of the 
primary system. By following the stability analysis presented in 
section 3.2.2, this small amplitude solution is found to be 
stable. The existence of such a stable, small amplitude solution 
is also confirmed by obtaining it through direct numerical 
integration. It was also observed that the other impacting 
solutions, different from this small amplitude solution, obtained 
from equations (3.14), (3.7) and numerical integration match very 
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closely. 

Selecting the design frequency as , the parameters r^ 

Ir 

and d^ are chosen satisfying the relation 

r d = 1/ . (3.19) 

mo P 

The frequency responses obtained from equation (3.7) for a fixed 

value of r^ with d^ given by equation (3.19) are shown in Figures 

3.3(a) - 3.3(c) for various value of (j^. With the impact damper in 

operation, only the stable solutions are shown. It is seen that, a 

stable, small amplitude solution is always present around the peak 

frequency n irrespective of the value of . Apart from this low 

amplitude solution, a high-amplitude stable solution is also 

present. However, there is always a region around where only 

the low amplitude solution exists. We define the frequency-range 

around Q , where the impact damper introduces a single, low 
P 

amplitude solution in place of the resonance branch of the primary 

system, as the suppression band (AQ^). These suppression bands, 

indicated in Figures 3.3(a) - 3.3(c), widen with increasing 

(i.e., decreasing relative non-linearity). However, within the 

suppression band, the percentage reduction in the amplitude is 

relatively insensitive to the degree of non-linearity. A decrease 

in (j. reduces the value of n_. Consequently, for a given value of 
I P 

r^, the optimum value of d^ obviously increases with decreasing 
(see equation (3.19)) as can also be noted from Figures 3.3(a) - 

3.3(c). It is known from previous works [77] on linear systems, 
that near the resonance frequency, the loose mass undergoes two 
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1.2 1.9 2.6 
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Figure 3.3 Frequency response of primary system with and without 

the impact damper. stable branch (harmonic 

balance) with damper, . stable branch without 

damper, unstable branch without damper, h = 0.215 

,rl = 2000 ,r = 0.12. a) 1.5, d^ =0.69, b) u, = 

m to t 

1, d^= 1.2 , c) = 0 ,d^= 2. 
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symmetric collisions per cycle. To examine what happens in the 
presence of a non-linear restoring force, the phase diagrams of 
the primary system were also studied. A representative diagram is 
shown in Figure 3 . 4 which clearly reveals the existence of two 
symmetric impacts/cycle. 

From here until section 3.5, the primary system will be 
considered to have only a cubic restoring force, i.e., = 0. The 
frequency responses of the system obtained from equation (3.7) for 
different combinations of r^ and d^ satisfying equation (3.19) are 
shown in Figures 3.5(a) - 3.5(c). The results obtained by the 
direct numerical integration are also shown in the same figures. 
The small amplitude motion is vindicated by the direct numerical 
integration (revealing two symmetric impacts per cycle) . However, 
the high-amplitude solution, obtained by integrating the equation 
of motion, revealed more than two impacts per cycle. The 
frequency regions having different number of impacts per cycle 
were also seen to be separated by narrow chaotic regimes. 

The number of impacts per cycle is not an important issue in 
the method of harmonic balance. For a single term approximation to 
yield meaningful results, only the energy content in the 
fundamental frequency is of primary importance. Of course the 
stability analysis carried out in this work is approximate and 
therefore, is not able to capture all the delicate instability 
regions. However, both the amplitude and the width of the 
suppression band are estimated accurately by the harmonic balance 
method which goes a long way to provide an initial design. Both 
the harmonic balance and numerically integrated results showed 
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n 

Figure 3.5 Frequency response of the primary system with and 

without the impact damper. stable branch (harmonic 

balance) , o-numerical integration, — •— stable branch 

without damper, unstable branch without damper, 

A - one-third subharmonic without damper, x one-third 

subharmonic with damper, h = 0.215, 0, rl = 

2000.(a) r = 0.12, d = 2; (b) r_= 0.25, d = 1 ; (c) 

m o m o m 

= 0.4, d = 0.625. 
o 
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that the motions of the primary and the loose mass are almost in 

anti-phase for and almost in phase for n>n . 

P P 

Let us define the set of values of the pair (r , d ) 
satisfying equation (3.19) as the optimum set. In the following 
discussion we vary (r^,d^) only in that set. A comparison of 
Figure 3.5(a) with Figure 3.5(b) reveals that the suppression band 

widens with increasing values of r . Thus, a high value of r is 

in m 

to be preferred for a wider suppression band. Moreover, higher the 

value of r^, lower is the value of the oscillation amplitude in 

the suppression band. But too high a value of r^ makes the system 

more prone to subharmonic oscillations (within the suppression 

band) as indicated in Figure 3.5(c). 

To determine the optimum values of r^ and d^, we compute 

the response of the primary system at the design frequency for 

three different values of r_ with varying values of d^. The 

maximum displacement of the primary mass is obtained both from 

equations (3.7) and numerical integration. These results are shown 

in Figure 3.6. The numerically integrated results are seen to 

yield a little higher value of the optimum d^ as compared to that 

obtained from the harmonic balance method. However, the motion of 

the primary mass at is seen to be rather insensitive to the 

value of r if the optimum value of d is used. Thus, depending on 
in o 

the required suppression bandwidth and available space, a value of 

r can be chosen and the first estimate of d_ should be obtained 
m o 

from equation (3.19). Around these values of r^ and d^, the exact 
optimum values can then be obtained from the numerically 
integrated results. The insensitivity of the primary response near 
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the optimum design ensures satisfactory performance of the impact 
damper, even if the actual parameters are a little away from the 
theoretical optimum values. These observations remain valid for 
different values of the parameter h as can be seen from Figure 
3.7. Furthermore, one can conclude that the HBM overestimates the 
optimum performance. 

Another important feature of . the impact damper for a 
non-linear primary system is the suppression of the one-third 
subharmonic response. Figure 3.2 reveals that the primary system 
without the damper can have one-third subharmonic oscillation near 
the design frequency. Figures 3.5(a) and 3.5(b) clearly show that 
the one-third subharmonic motion can even be completely suppressed 
by a properly designed impact damper. With higher values of r^^, 
the one-third subharmonic motion may appear but with considerably 
lower amplitude. Such control of the one-third subharmonic 
response is prominently displayed with a low value of h as 
indicated in Figure 3.3. 


3.4 Inelastic Impact Damper 

A mathematical model of an impact damper undergoing 
inelastic collisions is shown in Figure 3.9, where the impact is 
modelled by a linear spring with stiffness K 2 and a dashpot with 
viscous damping coefficient C^- The non-dimensionalized equations 
of motion can be written as 



Figure 3.7 Comparative study of analytical and numerical optimum 

design for various values of the parameter h. h = 

0.107, Q = Q (= 2.8), r^= 0.18, - . ~ . -analytical, # 

numerical; r = .12, .... analytical, e numerical, h = 
m 

0.15, Q = Q (= 2.37), r = 0.18, - - - - analytical, g 

pm 

numerical; r = 0.12, analytical, A numerical. 
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Figure 3.8 Frequency response of the primary system with and 

without the impact damper, h = 0.107, r, = 2000, r = 

1 m 

0.12. stable solution without damper, - - - 

unstable solution without damper, - ♦ - analytical 
stable solution with damper, o numerical low amplitude 
solution with damper, • numerical high amplitude 
solution with damper, A one-third subharmonic solution 
without damper, e one-third subharmonic solution with 
damper. 
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+ h + zj + = cos nx (3.20) 

• • • • • 

and + ^2 ^^ 2 ' ^2^ "" (3.21) 

with 

*l(Z3,Z2) = r^^CZj) + 25^ M(Z2) Zj, (3.22) 

«2'"2'"2> = <V^n> *<"2> + »("2> ^2 <=•=='> 


and M (Z^) = U (22-^^) + U(-Z 2 -d^) . (3.24) 

As shown in the Appendix B, the damping ratio ^ is obtained in 

w 

terms of the coefficient of restitution (x) as 




ln(y) 

(l+rJCTi^ + (Inix))^] 


(3.25) 


ICo study the effects of the coefficient of restitution 
on the performance of the impact damper, equations (3.20) and 
(3.21) are numerically integrated to obtain the primary response. 
The frequency response of the primary system is plotted in Figure 
3.10 for X = 0.75. We observe from this figure that the 
suppression band with < 1 is wider than that with x = !• 
However, the response amplitude with inelastic collisions is 
marginally higher than that with elastic collisions. In general, 
higher the impact damping (i.e. lower the value of x) wider is the 
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Figure 3.10 Frequency response of an inelastic impact damper, x = 

0.75, r = 0.25, d =1, h = 0.215.- • - • - stable 
m o 

branch without damper, unstable branch without 

damper, stable branch with elastic damper 

[analytical], o numerical with elastic damper, A 
one-third subharmonic without damper, + numerical 
with inelastic damper. 
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suppression band, but larger is the response amplitude (in the 
suppression band) . A similar behaviour was also seen for a linear 
primary system [69]. Thus, a high value of Xr signifying almost 
elastic collisions, is preferable for optimum performance in case 
of a synchronous operation. 

The variations of the amplitude of the primary system 
with d for various values of r and x are shown in Figure 3.11. 
From this figure, it may be concluded that higher the impact 
damping, lower is the required value of d^ for optimum 
performance. Moreover, as the impact damping increases, the 
primary response becomes more insensitive to the variation in d^ 
(for a near optimal design) . 


3.5 Improvement of the Design Using the EPL Technique 
In the previous sections the harmonic balance method (HBM) 
was used for designing and analysing the performance of an impact 
damper. In case of a perfectly elastic damper, a simple 

approximate relation for the optimum design was established. 
However, this optimum design had to be improved upon using 
numerical simulations. Further, for an inelastic damper only the 
numerical integration could produce meaningful results. The 
results obtained using the HBM were found to be quite inaccurate 
and are not presented here. Considering the difficulty and the 
amount of time involved in numerical simulations, in this section, 
we make use of the powerful EPL analysis (discussed in the 
previous chapter) for obtaining a good initial design. It has been 
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mentioned earlier that the optimum performance of an impact damper 
corresponds to a periodic motion having two symmetric impacts per 
cycle. In the HBM analysis one really does not require any 
information regarding the number of impacts. However, in the EPL 
analysis a periodic solution is constructed considering the number 
of impacts per cycle (two for the present problem) . It is also to 
be mentioned that for the present analysis, a spring-dashpot model 
of an impact is not really necessary. Rather, the concept of 

coefficient of restitution is used in the subsequent analysis to 
relate the velocities before and after an impact. 

3.5.1 Periodic Solution 

The equations of motion of the oscillator between two 
consecutive impacts can be written in the following 

non-dimensional form 

* * ’ 2 3 

Z^+ h Z^+ Z^+ Z^ = cos (nx) , (3.26a) 

Zp = constant (3.26b) 

and the impacts occur at Z,- Z = ± d 

1 p o 

where Z_ = x_ / x 

p 2 ' o 

For analysing the motion with two symmetric impacts per cycle, we 

set X = 0 at the instant of an impact occurring with the right 

side of the container when the next impact occurs at x = (n/fl) 
with the left side of the container. In the subsequent analysis a 
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+ sign denotes the instant just after an impact and a - sign 
denotes the instant just before an impact. We concentrate only on 
the motion occurring between < r < By the equivalent 
linearization procedure, equations (3.26) are rewritten as 


• # • ry 

Z^- 2p Z^+ Z^+ Z = cos (Qt) (3.27a) 

Zp = -V ^ < T < ( 2 )■ (3.27b) 


where p = - h/2 and v is the constant velocity of the loose mass 
between two successive impacts. The solution of equation (3.27a) 
is written as 


Z^(T) = e^^ (Bj^ cos 7]x + sin tjt) + M cos (Or + (p + 9) 


4- TT - 

^ O < T < (g) 


(3.28) 


where 


^ - P / 


tan 9 


2 p n 



and 


1 

M = - 

/ (w^ - n^)^ + (2pn)^ 



with 
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The symmetry of the motion implies 

Z^(0) = - Z^( 2 ) = Z° (say) 

Z^(0'^) = - Zi( n (say) (3.29) 

and ^ ^ 2° (say) 

The conservation of momentum during the impact and the definition- 
of the coefficient of restitution imply : 

X = (V + Z^)/ (V - Z°) (3.30) 

2 d^ + 2 Z° = ^ (3.31) 

o u 

and 2 r„ V = Z^ - Z° (3.32) 

in 

where x is the coefficient of restitution and r^^^ = 

mass ratio. Using the above conditions we eliminate other 

constants to obtain 

f3 ( Ke, V ) = 0 (3.33) 

To determine we assume the secondary solution as 

Z* (T) = cos (Qt + 0^) 
when from equation (2.12) we get 


K 


e 



(3.34) 
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and the condition (2.13) gives 


A 


2 

1 


' 2n 

7T 

V 



(■c) 


O 


i^T 

e 


dT 


2 


(3 . 35) 


Thus, using equations (3.33) - (3.35) both the primary (Z^(t)) and 
secondary (Z^(t)) solutions for the displacement of the main mass 
are obtained. We may recall (mentioned in section 2.4) that the 
EPL provides the time response Z^(t) (equation (3.28)) accurately 
whereas the single term HBM given by equation (3.5) may fail. An 
example substantiating this statement, in the context of the 
system under' consideration, is elaborated in the Appendix C. 

3.5.2 Stability Analysis 

The stability analysis of the solutions discussed in 

Section 3.5.1 is carried out using the method of error propagation 

described in Section 2.2.2. However, the system under discussion, 

in contrast to those discussed in Sections 2 . 2 and 2.3, has two 

degrees-of-freedom. Consequently the required perturbation vector 

whose change over a cycle is correlated by the matrix [P] (see 

equation (2.24)) has four components and the matrix [P] becomes a 

4x4 matrix. The four components of the perturbation vector are 

* 

taken as (AZ^, A2^, Av, A^) where A\^ is the perturbation in the 
forcing phase. As usual, if all the eigenvalues of [P] lie within 
the unit circle around the origin of the complex plane then the 
solution is stable. 
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3.5.3 Results and Discussions 

The optimum values of obtained for a given design 

frequency via different methods are plotted in Figures 3.12(a) and 
3.12(b). The Stable regions for the motion (with symmetric, two 
impacts per cycle) are also delineated in the same figures. The 
following conclusions can be drawn from Figures 3.12(a) and 
3.12(b): 

(i) The optimal values of d^ obtained using the EPL technique are 
in perfect agreement with the numerical results. 

(ii) For a perfectly elastic damper, the optimum values may be 
unstable for some choices of r^ as seen from Figure 3.12(a). 

(iii) As shown in Figure 3.12(b), a slight impact damping 

stabilizes the motion corresponding to the optimum design in the 

entire range of the r_ considered in this work. 

in 

(iv) A relation like 

r^ d_ = constant 
m 0 

is valid only for an elastic damper. However, the value of the 
constant is different from 1/ which was obtained using the HBM. 
For an inelastic damper such a relation is not uniformly valid, 

the discrepancy being more pronounced for low values of r . 

m 

3.6 Overall Suggestions for the Design of an 
Elastic Impact Damper 

So far we have presented three different methods, viz., HBM, 
EPL and direct numerical integration for arriving at the design of 
an impact damper for controlling the high-amplitude vibration of a 
harmonically excited Buffing's oscillator. Each of these methods 




with two impacts per cycle. 




78 


has certain advantages and disadvantages as discussed below. 

The HBM not only provides a good initial design of an elastic 
impact damper with minimum effort but also gives the overall 
frequency response without bringing in the actual number of 
impacts occurring per cycle. However, The initial design needs 
refinement to arrive at the optimal design. The EPL technique 
provides the most accurate optimal values. This method is also 
capable of tackling an inelastic damper. But this method requires 
prior knowledge of the actual number of impacts per cycle. With 
increasing number of impacts, the algebra involved in this method 
becomes unwieldy. Numerical integration being costly in terms of 
computation time should be done to check the final design. 

Considering the merits and demerits of various methods, it 

can be suggested that an initial design should be made (for a 

given Qp) following equation (3.19). First the optimum value of d^ 

should be obtained with a low value of r . This optimum value 

m 

should then be refined using the EPL technique. The frequency 
response of the system with these optimum values should be 
obtained using the HBM. If the suppression band needs to be 
widened, a higher r^ has to be tried. Once the final design is 
decided, direct numerical integration of the equations of motion 
has to be carried out to ensure the non-existence of any 
high-amplitude solution (not covered by the HBM, e.g., 
subharmonics.) in the frequency range of interest. 
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3.7 Use of an Impact Damper for Guided Branch Transition 

The concept of a guided branch transition, or migration 
control as it is called, is recent in the literature [103,104]. 
Non-linear systems generally have multiple, steady-state, periodic 
or strange (chaotic) attractors for fixed parameter values and all 
such attractors coexist with their respective basins of 
attraction. The role of a migration controller is to lead the 
system to a pre- assigned attractor from any initial condition. If 
the initial condition falls within the basin of an undesirable 
attractor, then the migration controller operates and guides the 
system to the basin of attraction of the desired attractor. Once 
this is achieved, the controller is released. The whole concept 
can be recast mathematically as follows: 

A forced non-linear oscillator can be described by an autonomous 
first order vector differential equation as follows: 

Z = f (Z, p(e)) 

8 = a (3.36) 

where Z e IR and 0 e ; 

f : X , 

1 I 

p : S -» R with I s n , 


the dot represents differentiation with respect to time x and a is 

a constant. With control variables Z„ and the control force 

c 

f (Z ) , the system equations take the following form: 

c o 



sa 


z = f (z, p(e)) + , 

Zc= f* (9, Z, Z^) 


(3.37) 


and 

e = 

a 

where 

Z^ e 

c 



^c = 

r’^ -> R^ 

and 

* 

^c ^ 

r” X r”* 




We represent the basin of attraction of the undesirable attractor 

A , of equations (3.36) by the set B. and that of the attractor 

ud 

of equations (3.37) by . Let the basin of attraction of the 


goal attractor A, of equations (3.36) be represented by B. . Thus, 

a 


B^ £ X 

^Ud 


B, 


£ r” X 


and 




We define the flow represented by equations (3.36) and (3.37) by 


and r^, respectively. Let the initial conditions be 

represented by a point (Z^,8^) € B^ such that 

ud 

Urn rJ^(Pi>^(z^,eo)) ^ Au^. 

r 00 

But if the control force is added, gets extended to 

P, (Z , e , Z ) € B- and 
1' o' o' co' A 

c 

lim r^(Pi(z^, z^^)) ^ A^£ B^ . 

T — > 00 C 

If the control is released when the motion reaches a point P 2 € 

[A^-n B^ ], then P^ is restricted to P 2 (Zj^, 9^^) [e B^ ] when 
^d d 

®l)) -- B ]. 

T — > 00 d 

The whole control action is schematically shown in Figure 3.13. 
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In Sections 3.2 and 3.4, we have seen that an impact 
damper can reduce the response of a harmonically driven, Duffing's 
oscillator, especially in a frequency-band corresponding to the 
high-amplitude resonance branch. Using the multiplicity features 
of the steady-state motion of the Buffing's oscillator, here we 
use an elastic impact damper as a migration controller. It will be 
seen how starting from the basin of attraction of the resonance 
solution, the system goes to the non-resonance branch when an 
impact damper is used as a controller. It is generally seen that 
around the higher end of the resonance branch, the probability of 
finding an initial condition for the non-resonance motion is 
higher near the trivial initial conditions. Around these trivial 
(zero) initial conditions, a wide zone exists which forms the 
subset of the basin of attraction of the non-resonance motion 
[105]. It has already been seen in sections 3.2 and 3.3 that an 
impact damper brings the primary system to a near-trivial motion. 
Thus, the points on the steady-state attractor of the system 
(including the impact damper) might be completely embedded within 
the basin of attraction of the non-resonance motion of the primary 
system (without the impact damper) . At , the steady-state, the 
primary motion might choose the initial point suitable for the 
non-resonance oscillation when the impact damper can be switched 
off. However, as the system is non-autonomous , the phase of the 
excitation at the time of withdrawal of the control is very 
important. The phase should correspond to the desired goal. Thus, 
using the migration control one can avoid the continuous noisy 
operation of an impact damper. 



S3 


A numerical example showing the feasibility of the above 
scheme is given below: 


Say 



and 


f 




a = Q 


where f^ (Z^,Z^,a) 



cos & 


If an impact damper is used as a controller, then 


Z 


c 




^c 




and f%(e, Z, Z^) = ^ 1_, 

where ‘^o^ ^ ^^ 2 " *^ 0 ^ 


+ (Z 2 + d^) U (-Z 2 - d^)|. 


Let us take r^ = 2000, h = 0.215, = 2 and r^ = 0.12. 

We consider the operation of the impact damper as a migration 
controller at three different values of the excitation frequency. 
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(i) n = 1.95 

with the above mentioned functional form and prescribed parameter 
values, we numerically integrate equations (3.36) staring from an 
initial condition. Let the initial condition be within the basin 
of the resonance motion. After the motion has reached the steady- 
state, we make the controller (i.e.,the impact damper) on and 
again numerically integrate equations (3.37) up to the 
steady-state. When the control is switched off at some point, the 
subsequent motion leads to the non-resonance branch. We repeat the 
procedure to ensure that even if the control is released at any 
arbitrary forcing phase, the motion still migrates to the 
non-resonance branch. One example showing the primary response 
(Z^) plotted against the non-dimensional time (x) is shown in 
Figure 3.14. The instants of switching on and off the impact 
damper are indicated in this figure. 

(ii) n = 1.7 

The steady-state attractor of the system with the controller 
(equations (3.37)) passes through the basin of one-periodic 
non-resonance motion and three-periodic motion of the system, 
described by equations (3.36), as shown in Figure 3.15. So 
depending on the time of release of the control, the system might 
choose initial conditions (definitely on the above mentioned 
attractor) leading to either of the two states. Thus, one may need 
several "on-off's to reach the non-resonance one-periodic motion 
as shown in Figure 3.16. 

(iii) £1 =1.85 

The response is plotted against time in Figure 3 . 17 where we 





Figure 3.14 Oscillatory signature of the primary system with 
migration control, n = 1.95, h = 0.215, r = 0.12, 
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Figure 3.15 Steady-state attractor of the primary system. SI = 1.7, 

r = 0.12, d = 2, h = 0.215. A represent the points 
m o 

from which the motion without the impact damper goes 
to one-third subharmonic oscillation. X represent the 
points f rom which the motion without the impact damper 
goes to the non-resonance motion. 





Figure 3,16 Oscillatory signature of the primary with migration 
control. Q = 1.7, other parameters are same as in 
Figure 3.14. 
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release the control before the steady state is reached. The 
motion is seen to reach a three-periodic solution. Several 
switchings of the controller ultimately lead the system to the 
non-resonance one- periodic motion. 

For a number of parameter values tested, it was observed 
that several on-off switchings spanning, in total, not more than 
350 cycles of operation could lead the response to the 
non-resonance branch. For a safe application of an impact damper 
as a migration controller, one needs to design the impact damper 
for its optimum performance. Thereafter, a hit and try control 
algorithm depicted in Figure 3.18 can be followed for migrating to 
the desired motion. 


3.8 Conclusions 

The performance of an impact damper attached to a 
harmonically excited Duffing' s oscillator has been presented. The 
existence and stability of a low amplitude motion, near the 
high-frequency, large amplitude resonance branch of the primary 
system, has been established using the harmonic balance method and 
the EPL technique. The optimum design predicted by the harmonic 
balance solution is refined by the EPL technique and numerical 
integration. For a given mass ratio (r^) , the optimum free space 
for the loose mass (d^) is obtained. The primary response is seen 
to be rather insensitive to the values of both r^ and d^ near this 
optimal design. With the optimum design, two symmetric collisions 
per cycle are always seen to occur. The suppression band of the 
impact damper narrows down with increasing non-linearity in the 
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restoring force, but the percentage reduction of amplitude at the 
design frequency is seen to be insensitive to the degree of 
non-linearity. An impact damper is very effective for reducing (or 
even eliminating) the one-third subharmonic response. An inelastic 
damper widens the suppression band at the cost of higher primary 
response. An "on-off” migration control strategy, using an impact 
danpeJ^/ ^ guided transition from one solution branch to the 
other is proposed. The feasibility of the proposed scheme is 
verified through numerical examples. This scheme, applicable to 
other non-linear primary systems having multistable solutions, 
obviates the need of continuous, noisy operation of an impact 
damper. 



CHAPTER-4 

IMPACT DAMPERS FOR CONTROLLING SELF-EXCITED OSCILLATION 

4.1 Introduction 

If the vibration of a system is excited by a force which 
depends on the state variables representing the motion, then it is 
called self-excited vibration. Self-excitation is a common source 
of vibration in various physical systems, such as (galloping of) 
transmission lines, (chatter in) metal cutting machine tools, a 
rotor bearing combination and a number of dry friction driven 
oscillators. The mechanism of self-excitation is identified with 
the presence of a non-linear function of state-variables in the 
mathematical model of the system. The most widely used foma of 
this function is van der Pol type damping which has a negative 
linear part destabilizing the equilibrium position of the system 
and a non-linear part limiting the growth of the instability. As a 
result, a new dynamic equilibrium (the so called limit cycle) is 
approached and the system continues to vibrate autonomously. Of 
course, a self-excited system can also be simultaneously subjected 
to some other external forces when the system is termed 
non-autonomous or forced. 

Several attempts have been made to quench the self-excited 
vibrations in physical systems. Mansour[106] has investigated the 
role of a conventional dynamic vibration absorber (DVA) in 
suppressing the self-excited oscillations. He investigated the 
possibility of using a linear DVA for stabilizing the 
self-excitation of a van der Pol oscillator and found the results 


92 


93 


not very promising. Tondl [107] extended the analysis of Mansour 
for the same system and concluded that the absorber can be 
effective only if the mass ratio (ratio of the absorber mass to 
that of the main system) is very high. Tondl [108] has also 
studied the performance of a resilient foundation mounting towards 
suppressing the self-excited vibration of a rotor system. He 
investigated the role of dry friction [109] in controlling the 
self-excited vibration and found that dry friction may be 
effective in stabilizing the equilibrium position of a van der Pol 
type oscillator under soft perturbations. Recently, Asfar [110] 
has proposed the use of a Lanchester damper for quenching both the 
autonomous and forced self -excited vibrations of a single 
degree-of -freedom system. He concluded that the tuning of the 
damper is not important for its efficient operation. 

The major emphasis in all the above mentioned works was on 
the stability of the equilibrium. For successful design of these 
devices, an accurate knowledge of the strength of the 
self-exciting force was necessary so that a proper amount of 
external damping can be incorporated to compensate the negative 
damping (which is the source of instability) . In this chapter we 
investigate the performance of an impact damper for controlling 
both free and forced vibrations of a self-excited oscillator. 
Complete quenching (equivalent to the stabilization of the 
equilibrium) is obviously not possible by an impact damper, since 
an impact damper is operative only when there is oscillation. A 
van der Pol oscillator is analysed and numerical results are 
presented to show the effects of various parameters. 
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4.2 Equations of Motion 

A mathematical model of a single degree— of— freedom van der 
Pol oscillator with an impact damper is shown in Figure 4.1(a). 
The differential equations describing the motion of the system 
between two consecutive impacts can be written in the following 
non-dimensional form 



+ e 



1 ) + 


Fq cos(Q t) 


(4.1a) 


• • 

and 2=0 (4.1b) 

Jr 

^ I - 2p I < do 

with Wq =(K^/m^)^/^, ■„ Xq = v'g^, Q = (o/w^ , F^ = F/(m^XQWQ) 
e = e^XQ / (m^WQ^ ' ^1 "" ^p ^ ^2^ ^0' d/2XQ, t = w^t 


and the dot denotes differentiation with respect to r. 

When I - Zp | = a collision occurs satisfying the following 

conditions : 


Zi = «i + 

1 Ip 


* + * - 

Z. = a- Z^ + 
1 3 p 


a. 



a. Z^ 
4 p 


where the superscripts + and - denote, respectively, the 
conditions after and before impacts. Using the conservation of 
total linear momentum and the definition of the coefficient of 




(b) 


Figure 4.1 van der Pol oscillator with an impact damper, (a) zero 
contact time, (b) finite contact time. 
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restitution it is easy to obtain s as follows: 


“1 

= {X 

- r )/ 

( 1 

+ 

X) 


“2 

= (1 

m' 

(1 

+ 

X) 


“3 

= X ( 

1 + r 

la 

)/ 

( 

1 -1- 

X) 

“4 

= (1 

• r y 
m 

)/ 

( 

1 + 

X) 


where x is the coefficient of restitution and r = m^/ m, , is the 

m 2 1 

mass ratio. 

An alternative way of modelling an impact damper would be 
through the introduction of a spring and a dashpot as shown in 
Figure 4.1(b). The differential equations describing the .motion of 
such a system are written in the following non-dimensional form: 

h ^ ^^1 ^ + ""l ^ 1 (^ 2 ) + 

and (4.2a) 

z; + r,*,(Z2)/r^ + *2<"2)22 ' (4.2b) 

where 

*1(2^) = (Zj - do)U(Z2 - dj,) + (Zj 4 d„)U(-Z2- dj,) 

* 2 (^ 2 ) = “<22 - ‘»0 ) + "<-22 - > 

when U(.) is Heaviside's step function defined as 

U(C) =1 ^ C > 0 

= 0 -v- c < 0 

and r^ = K^/ = (x^ - -x .^) and the damping ratio is 

obtained in terms of the coefficient of restitution x used in the 
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first model as (see Appendix B) 


C = C- 
2 




ln(x) 

{(1 + r^) ( TT^ + (ln(;t))^)} 


(4.3) 


This expression is valid except in a small neighbourhood of x = 0. 

With a very high numerical value of r^ (~ o(10^)) in 
equations (4.2), the results obtained from both the models 
described above are found to be quite close. However, it has been 
seen that the second model described by equations (4.2) offers 
special advantages over equations (4.1) during numerical 
integration. 

It may be mentioned here that an impact damper (shown in 
Figure 4.1(b)) is a special form of a non-linear dynamic vibration 
absorber (DVA) where the non-linearity is introduced through the 
clearance. It can be easily seen that the equations of motion 
reduce to those of a linear DVA if d^ in equations (4.2) is set to 
zero. 


4.3 Autonomous vibration 

For the analysis of autonomous vibration, we put F^ = 0 in 
equations (4.1) and (4.2). It is a well known fact that an impact 
damper, used to suppress the vibration of forced linear/non-linear 
oscillators [77,88], operates efficiently with two symmetric 
impacts occurring in a forcing cycle. Assuming the same to be true 
during one cycle of a self-excited oscillator, we shall construct 
a periodic solution with two symmetric impacts per cycle and 
analyse its stability. 
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4,3.1 Periodic Solution 

In this section we present a syminetric periodic solution with 
two impacts per cycle. One half of such a cycle is depicted in 
Figure 4.2. Without any loss of generality, let us arbitrarily set 
the time origin (x = 0) at the instant of an impact with the right 
hand side wall of the container. Now from the symmetry condition, 
it is known that the next impact will occur at x = where 0^ 
is the unknown frequency of the limit cycle. First we transform 
equations (4.1) into equivalent linear forms given by 


+ 2 p = 0 


{4.4a) 


and 


Z = - V 
P 


(4.4b) 




- "p 


< dg and 0^< x < (tt/Q^) 


where p is an unknown constant to be determined later and v is the 
velocity of the loose mass. From the symmetry of the motion we 
have the following conditions; 

Zi(0) = - Z^(7r/np = Z° (say) 


Zj^(0‘^) = - Z^(Tr/n^)'^ = Z^ (say) 


and Z 


10 


^(TT/n^) = - Z^{0 ) = - Z (say) 


(4.5) 


The definition of the coefficient of restitution, equation (4.4b) 
and the conservation of total linear momentum yield 
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n/Q^ 


Figure 4.2 Half-cycle of a typical periodic solution having two 
symmetric impacts per cycle (solid line) and its 
arbitrarily perturbed form (dashed line). 
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;t=(v + z^)/(v-z°) , 

2 dg + 2 Z° = V 7T /n^ 

* 1 * 0 

and 2rv=Z-Z. (4.6) 

m ' ' 

Solving equation (4.4a) one can write the motion between two 
consecutive impacts as 

Z^(T) = exp(-p T) ( d^ cos ( 7} r ) + d^ sin ( tj t) ) (4.7) 

where tj = 1 - 

Now using conditions (4.5) and (4.6) one can eliminate other 
. ■ * 1 * 0 

unknowns viz., d^^, d 2 , v, Z , Z to get an equation of the form 

f^( p , ) = 0 (4.8) 

For determining p let us consider a single term Fourier series 
approximation of Z^(t) given by 

Z*(t) = a cos (fi^T) (4.9) 

satisfying the equation 
27T/n, 

r ^ 2 2 

j C Z* ( Z* - 1) - 2 p Z* J- dT = 0 (4.10) 

0 

Following this procedure p turns out to be 


p = c ( a2/4 -1 )/2 


(4.11) 
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It is to be noted that equations (4.7) and (4.9) are two different 
representations of the same solution. Thus, these two should be 
equal in some least square sense which implies that the first 
Fourier component of (4.7) should be given by (4.9) i.e., 


2n/Q 


n 


(T) 


(4.12) 


It may be mentioned that d^ and d^ can be written as explicit 
functions of the unknown a. Thus, Z^(r) is also an' explicit 
function of a and n^. So equation (4.12) is finally transformed to 
the form 


a, ) = 0 (4.13) 

Now using equations (4.8) ,(4.11) and (4.13) one can find out a 
and when the other unknowns can also be easily determined. In 
fact, periodic solutions having any number of impacts can be 
constructed in a similar fashion. 

4.3.2 Stability Analysis 

For analysing the stability of the solution obtained in the 
previous section, the method of error propagation explained in 
reference [69] has been used. A brief sketch of the method is 
given below. Referring to Figure 4.2, let the solid line represent 
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the half-cycle of a constructed solution whereas the dotted line 
represents its perturbed form. The initial perturbation vector 
is represented by 

T 

= I AzJ; , Lzl , Av^, At^ I 

This perturbation vector at the end of a half-cycle is given by 

T 

Cf = j AZ^ , AzJ , Av^, At^ J- 

A linearized relation between these two vectors can be written in 
the closed form as 

Cf = C P ] Ci (4-14) 

where [P] is a 4x4 matrix whose elements are functions of the 
system parameters. As the solution is symmetric, the propagation 
rule for the perturbation vector during the next half of the cycle 
is also given by equation (4.14). Thus, the stability of the 
solution will be governed by the eigenvalues of the matrix [P]. If 
all the eigenvalues of the matrix [P] lie within the unit circle 
around the origin of the complex plane, then the solution is 


stable. 
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4.3.3 Results and Discussions 

In this section, the role and design of an impact damper for 
reducing the limit cycle amplitude of a van der Pol oscillator is 
reported. Before going into the detailed discussions, one should 
remember that the van der Pol oscillator has a non-dimensional 
limit cycle amplitude 2 irrespective of the value of e {<1) . So 
our aim is to reduce this limit cycle amplitude by properly 
designing the impact damper which is described by the parameters 
r^jj, d^ and It is obvious that a lower value of x implying 

larger amount of dissipation counteracts the negative linear 
damping and thus helps to reduce the amplitude of the oscillator. 
Therefore, a very low value of x (= 0.2) is chosen. A material 
pair having so low a coefficient of restitution is not very 
uncommon in practice [77]. To determine the effects of the other 
parameters, viz., d^ and r^, results are obtained both 
analytically and by direct numerical integration. The integration 
is carried out using fourth order Runge-Kutta-Merson algorithm 
with adaptive step size control. As mentioned earlier, integration 
of equations (4.1) is slightly difficult, since it involves the 
determination of the exact instant of an impact. However, 
equations (4.2) are easy to integrate due to the presence of a 
finite duration impact. In fact, the amplitudes of the limit 
cycles and their stability characteristics obtained from both the 
models (with a suitable choice of r^) are found to be quite close. 

The variation of the stable limit cycle amplitudes, 
with d^, corresponding to the symmetric motion with two impacts 
per cycle, is plotted in Figures 4.3(a) and 4.3(b). It is seen 
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that a significant reduction of the amplitude is possible with 
proper choices of r^^^ and d^. It may be noticed that for some 
values of r^, there is a range of d^ where two stable limit cycles 
(both having two symmetric impacts per cycle) may coexist giving 
rise to a so called 'jump' phenomenon. 

The results obtained from direct numerical integration are 
also plotted in Figures 4.4(a) and 4.4(b). It should be mentioned 
that, when the numerically integrated results match those 
predicted by the analysis, the motion is also seen to consist of 
two symmetric impacts per cycle. The extra results obtained by 
numerical integration reveal altogether different and complicated 
types of motion. However, the major objective of this work is to 
determine suitable values of r^^ and d^ for a good design of the 
impact damper. Therefore, we shall not discuss in details all 
these complicated motions which may be one of the following types; 

i) stable periodic solution with asymmetric two impacts per cycle 
and stable periodic solution with more ( than two) impacts per 
cycle, 

ii) quasiper iodic, 

iii) intermittent and chaotic. 

One may conclude from Figures 4.3(a) and 4.3(b) that, for a 
good design of the impact damper, r^ and d^ should be chosen so as 
to minimize the amplitude of the limit cycle. However, if for the 

design value of r , there exists a 'jump' region in the 1 Z, 1 vs 

d^ plane, then d^ should be chosen outside the jump region. This 
design procedure presumes an exact knowledge of the parameter e 
which may not be possible in a real life situation. In fact, the 
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value o^ 'l^lixnax ^ pair of (r^, d^) values is highly sensitive 

to the value of c as can be explained from Figure 4.4(a). For 

example with = 0.1 the optimum d^ based on c = 0.2 is around 

1.8. If instead, the actual value of e is 0.25 then at the same 

(r , d-) combination, the value of | Z. I increases by a factor 

of 2 . On the other hand if the actual value of e is 0.15, then the 

desired two symmetric impacts per cycle motion does not exist at 

all. The variations of |Z 1 with e for different (r ,d„) 

combinations are shown in Figure 4.5. Finally we conclude that the 

design procedure should not be restricted only towards minimizing 

|Z, I but should also take into account its sensitivity to the 
' I'max 

value of c . 


4.4 Forced Vibration 

In this section, the performance of an impact damper for 
controlling the vibration of a forced, van der Pol oscillator is 
discussed. A detailed account of the dynamics of a van der Pol 
oscillator can be found in the literature [ill]. It is well known 
that the phenomenon of frequency entrainment occurs if the 
external forcing amplitude is greater than a critical value and 
the forcing frequency is close to an integral or sub multiple of 
the frequency of the free, limit cycle oscillation. As a result, 
the oscillator resonates at the forcing frequency. All types of 
resonances viz, superharmonic, harmonic and subharmonic 
entrainments are theoretically possible. However, the harmonic 
entrainment region extends over a large frequency band-width and 
in this zone the oscillator vibrates with a large amplitude. 






Therefore, the discussion is limited to the control of the forced 
vibration in the harmonic entrainment zone. As before, we assume 
that the motion with two symmetric impacts per cycle offers the 
best mode of operation and proceed to construct the periodic 
solution corresponding to this mode. A stability analysis is also 
carried out. 

4.4.1 Periodic Solution 

For the construction of a symmetric motion with two impacts 
per cycle we refer back to Figure 4.2. The time origin r = 0 is 
set arbitrarily at the instant of an impact with the right hand 
side wall of the container. Consequently an unknown phase ^ is 
introduced in the forcing function. The motion of the oscillator 
between r = o"*" and x = (tt/Q) is described by the following 
equivalent linear differential equations 
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/ (1 - + (2pn)^ 

and e =; tan“^ ( 2pQ/(l - n^)). 
Conditions (4.5) and (4.6) 
conservation and the coefficient 
with Eliminating d^ , d^, 

the form 


expressing symmetry, momentum 

of restitution still hold good 

• 1 *0 

(p, Z , Z we get an equation of 


i V, p ) = 0 

As before, we consider 


to obtain 


Z*{r) 


a cos (Q T + 0^) 


(4.16) 


p = c ( a^/4 -1 )/2 (4.17) 

Moreover, equation (4.12) remains valid and we finally h: ve a 
relation like 


f 2 ( V , a ) = 0 (4.18) 

Now using equations (4.16) - (4.18) one can determine a and v. 

Thereafter, the other unknowns are obtained easily. 

The stability analysis of the solutions so constructed is 
carried out by the method of error propagation explained in 


section 4.3.2 
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4,4.2 Results and Discussions 

For the non— dimensional van der Pol oscillator under 
discussion, the limit cycle frequency is approximately unity. 
Thus, the fundamental entrainment takes place at and near n = 1. 
So we set our design frequency at n = 1. For most real life 

applications r^^^ < 1. It is well known, in the context of the 

impact damper for a forced passive oscillator, that the maximum 
attenuation is possible for perfectly elastic impacts (x = 1 ) • 
However, in the present case, the motion with two impacts per 
cycle is found to be unstable for perfectly elastic impacts. 
Therefore, we set x to be less than unity such that the motion 
with two impacts per cycle becomes stable. A parametric study 
showing the effects of the design variables, viz, x, d^ and r^^ on 
the amplitude of the forced vibration and the frequency band of 
suppression is reported. The results are obtained both 
analytically and by numerical integration. 

The variations of the amplitude of vibration, with d^ for 
various combinations of r^ and x , are shown in Figure 4.6 at n 
=1. At this frequency, the oscillator without the impact damper 
under the same level of excitation showed an amplitude greater 
than 3 which clearly demonstrates the usefulness of the impact 
damper. From Figure 4.6 it is obvious that an optimum d^ exists, 
the value of which depends on the values of r^ and x> The two 
impacts per cycle motion remains stable over a range of d^ around 
this optimum value. The graphs also suggest that just as in the 
case of an impact damper for a passive linear or non-linear 
oscillator [88], the value of x should be as high as possible. The 



Figure 4.6 Variation of the amplitude of forced vibration with 

d^. Fq = 0.5, n = 1.0, e = 0.1. Canalytical), 

A (numerical) for 0.25,;^; = 0.5), 

(analytical), o (numerical) for 0.25,;^; = 0.7). 

(analytical), n (numerical) for (r = 0.15,;:^ = 
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variation of optimum with for given values of the other 

parameters is shown in Figure 4:7. The limits on the parameter 
values maintaining the stability of the motion with two impacts 
per cycle are also indicated in the same figure by the dashed 
lines. The graph of Figure 4.7 fits the relation 

'^0>opt ' 

quite accurately. It has been found that at n = 1, the following 
approximate relation can be written 

where (3 ~ 0.06, X ~ 1.85, y - 1 , S ~ 0.12. 

The frequency response characteristics of the oscillator 

around Q =1, without and with an optimally designed impact damper 

are shown in Figure 4.8. The design of the impact damper may be 

optimized either at n =1 or at the operating frequency (around fJ 

=1) . In the latter case, of course, the value of d^ has to be 

actively controlled by sensing the excitation frequency. The 

frequency range over which the motion of the oscillator (with two 

impacts per cycle) remains stable is also indicated in Figure 4.8. 

This range of frequency (AQ ) is known as the suppression band. As 

s 

expected,, the actively controlled impact damper performs better 
(than the one passively tuned at £2 = 1 ) both in terms of 

amplitude reduction and suppression band. 

Numerical results for the passively tuned (at £2 = 1 ) impact 
damper suggested that higher the value of r^ (of course, 
satisfying equation (4.19)), not only the reduction in amplitude 
is more, the suppression band is also wider. On the other hand, if 





Figure 4.8 Frequency response of the oscillator. F^ = 0.5, x = 

0.7, e =0.1, r = 0.25, d. = 2.9. without 

m u 

impact damper. - ' • - ' * - with damper (passive, 
analytical), A (passive, numerical). 

(active, analytical), o (active, numerical). 




the coefficient of restitution is increased, then the reduction in 
amplitude is more at the cost of narrowing the suppression band. 
These facts are illustrated in Figure 4.9. The statements made so 
far regarding the performance of the optimally designed impact 
damper are true for given values of and e. However, the values 
of Fq and c in real life may not be known accurately. The 
sensitivities of the performance of an optimally designed impact 
damper (for synchronous operation at fj = 1) with the variation in 
Fq and e are indicated in Figure 4.10. It is observed that a 
higher value of r^ takes care of wider variations in the forcing 
whereas a lower value of x takes care of wider variations in the 
value of e . 


4.5 Conclusions 

For some combinations of the mass ratio (r ) and free space 

m 

(d^) parameters of the impact damper, an autonomous van der Pol 
oscillator can have more than one stable limit cycles giving rise 
to 'jump' phenomenon in the amplitude-dQ plane. For efficient 
operation, the impact damper should have a low coefficient of 
restitution and undergo two impacts per cycle. The parameters 
should be chosen outside the jump region. For the forced 
oscillator, a high value of the coefficient of restitution is 
desirable. An optimum design should attempt at both reducing the 
amplitude and widening the suppression band. An actively 
controlled damper performs better than a passively tuned one. The 
optimum parameters are sensitive to the levels of self and forced 
excitations. 


Figure 4.9 Optimum frequency response characteristics for various 

values of r and >;• 0.5, e = 0.1. 

m 

(analytical), o (numerical) for (r^= -25, x = 

0.5, = 3.1) . (analytical), d (numerical) for 

(r = 0.19, x=-'^ ~ * (analytical), A 

(numerical) for ( r = 0.25, x = ^*7, = 2.9). 
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4.10 Sensitivity of the designed performance to the 
variations of forcing amplitude and self-excitation 
parameter e. 0.25, c = 0.1, Q = 1, ~ 

0.7, dQ= 2.9 ) and ( r^= .19, e = 0.1, n = 1, 

= 0.7, d^ = 3.8) : sensitivity to the variation of F^. 
(r = 0.25, F = 0.5, Q = \,x = 0.7 d = 2.9) and 

m 0 u 

( r = 0.25, F = 0.5, Q =1, = 0.5. d = 2.9) : 

m O 

sensitivity to the variation of v. 



CHAPTER-5 

BIFURCATIONS AND CHAOS IN AUTONOMOUS SELF-EXCITED 
OSCILLATORS WITH IMPACT DAMPING 

5.1 Introduction 

Impact dampers in various forms are used to control the 
response of a vibratory system subjected to either forced or self 
excitations [70-72]. The response to a forced, harmonic excitation 
in the presence of an impact damper has been investigated by a 
nximber of researchers [67,77,86,91]. Masri [68,69] has constructed 
a general analytical solution using the piecewise linearity of the 
equation of motion. A method of analysing the asymptotic stability 
of the solution was also provided which was later refined by Bapat 
et al. [23]. A bifurcation study of the same system has recently 
been reported by Sung et al.[101]. Different kinds of local 
bifurcations e.g., pitchfork, saddle-node, Hopf etc., and global 
bifurcation such as homoclinic tangency have been studied. A 
period-doubling sequence leading to chaos has also been found. 

In a self-excited system, the exciting force is dependent on 
the motion of the system. The most common mathematical model is 
given by the autonomous van der Pol equation in which the position 
dependent damping force changes sign once in every cycle and acts 
as the source of energy for maintaining the vibration. The use of 
an impact damper for suppressing self-excited vibration is not 
uncommon. Both the transfer of momentum to the loose mass (from 
the main system ) and the dissipation of energy in every impact 
provide the mechanism for controlling the vibration. In this 
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chapter, we present a study on the persistent bifurcation 
structures and chaotic motions which remain invariant to the exact 
f\inctional form of the self-exciting force. Both analytical as 
well as numerical methods have been used. The analytical method is 
in line with that of Masri et al.[69]. The basic motivation behind 
undertaking this investigation is explained below. 

The first step, towards the understanding and prediction of 
the dynamical features of any physical system, is to construct a 
mathematical model. No real life system can be modelled in its 
entire complexity accounting for every detail. Moreover, 
idealizations are also necessary for obtaining a meaningful 
solution. However, certain phenomena predicted by such an 
idealized model may not be revealed under real life perturbations 
whereas some unpredicted phenomena may be exhibited by a real life 
system. Therefore, if possible, it is worthwhile to know which 
qualitative results of a mathematical model remain invariant under 
defined perturbations. The study of the qualitative invariance of 
the bifurcation structure is not very common in vibration 
modelling. This aspect has been considered in a recent work by 
Ravindra et al. [113]. It has been shown that the bifurcation set 
of a forced Buffing's oscillator remains topologically invariant 
under the possible perturbations of the index of the velocity 
dependent (i.e. dissipative), non-linear damping force. 

In the context of the present problem, it is the 
self-excitation mechanism which cannot be modelled accurately. No 
single model can prove to be the best and always represent the 
reality. Several variations of the model of the self-excitation 
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mechanism have been attempted. Of course, the self-excitation 
demands that certain topological conditions be satisfied by the 
function used to describe it. This topological requirement has 
been kept invariant in the analysis. The main aim is to detect 
those bifurcations which persist under defined perturbations, 
while some others vary with the model. 


5.2 Equations of Motion 

Let us consider a self-excited mechanical oscillator with a 
loose mass m 2 (the so called impact damper) as shown in Figure 
5.1(a). The collisions between the two masses, in most cases, can 
be adequately represented by identical, parallel combinations of 
spring and viscous damper as shown in Figure 5.1(b). This second 
model, allowing a finite time of contact, eventually smoothens the 


impact process and 

facilitates 

numerical 

integration 

of 

the 

governing equations 

of motion 

which can 

be written 

in 

the 


following non-dimensional form : 


+ e f(Z^, Z^) Z^ + Z^ + r^ 2 Cc ° 

(5.1a) 


Where the dot at the top indicates differentiation with respect to 
the non-dimensional time x (= w^t with t as the time) , 
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4 >i(Z 2 ) = (Z2- d^) U ( Z2 - dQ ) + (Z2 + dg ) U (■ 




^2(23) = U ( Z2 - dQ ) + U ( -Z2 - dQ ) 
and 

U(cr) = 0 ■, O' < 0 
= 1 , cr > 0 

with Z^= x^/ Xq, Z2 = (x^ - X2)/ Xq, Wq = 


?C = S/< 2 /^2 >' “^1 “ “o> - <*0 “ <2 =‘o>- 

In the above expressions, x^ with the dimension of length is 

suitably chosen depending on the form of the function f. The 

parameter C in equations (5.1) is given by [Appendix B] 
o 


In ( y ) 

y (1 + r^) ( It^ + (ln(;t:))^) 


(5.2) 


where x, the coefficient of restitution, is the ratio of the 
relative velocities of separation and approach during the impact. 
The self -exciting mechanism is inherent to the form of the 
fxinction t(Z^,Z^) , which has the topological configuration shown 

in Figure 5.2(a). A general expression of f, depending on a 
suitable choice for x^, can be written as 
either f(Zj^,Z^) = (a Z^ + Z^ - 1) , 


(5.3) 


or, f(Z^,Zj^) = (Z^ + ^ z\- 1) 

where a and 0 are bookkeeping parameters which enable several 
models of self-excitation to be handled by the form of f with a 
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and 0 as either zero or unity. The function f represented by 
Figures 5.2(b) and 5.2(c) can be represented by equation (5.3) 
with a = 0 and ^ = 0, respectively. The discontinuous function 
shown in Figure 5.2(d), also implying self-excitation, is amenable 
to exact analytical treatment. Moreover, this function also serves 
as a perturbation of the function shown in Figure 5.2(b) which 
represents the standard van der Pol form discussed so far in the 
thesis. 

A preliminary numerical simulation (reported in detail in 
section 5.5) of equations of motion (5.1) reveals the existence of 
several windows of periodic solutions as the value of the 
parameter d^ increases. By a periodic solution, we mean the 
solution with the minimum possible period which, in the present 
case, is approximately 6.6. Such a periodic solution is termed as 
the basic solution. Each window corresponds to a specific number 
of collisions per cycle. The number of impacts occurring in a 
period is seen to decrease with increasing d^. The simplest 
possible periodic motion observed during the numerical simulation 
consists ■ of two symmetric impacts per cycle. Other possible 
periodic solutions correspond to larger (such as 4,6 etc.) number 
of symmetric impacts per cycle. As the parameter d^ is varied, the 
periodic motion within each window undergoes a certain sequence of 
bifurcations and thereby gives rise to interwindow chaos. In order 
to determine the route to chaos, the bifurcation structure of the 
periodic solutions is investigated. In the sections to follow we 
construct these periodic solutions and analyse their stability. 
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5.3 Theoretical Analysis 

When the self-exciting function f is continuous, an 
approximate method based on the equivalent linearization can be 
employed to construct a basic periodic solution. For the sake of 
algebraic simplicity and relevance to the practical operating mode 
of an impact damper, the periodic solutions with two symmetric 
impacts per cycle are constructed, a detailed discussion of which 
with van der Pol damping has been included in chapter 4. The 
stability analysis of the solutions is also presented there. 
However, as stated earlier, an exact analytical treatment is 
possible for the piecewise constant form of f shown in Figure 
5.2(d). In what follows, this method will be briefly outlined for 
a syiametric periodic solution with four impacts per cycle. 

5.3.1 Periodic Solutions and Stability Analysis with Four Impacts 
per Cycle 

For the piecewise constant form of f, the periodic 

solutions with four symmetric collisions per cycle can be 

constructed as follows. One-half of a typical solution for is 

shown in Figure 5.3. In this figure, a point is defined as an 

impact node where an impact occurs. Similarly, a cross node is 

defined where a sudden change in the value of f takes place . Thus 

the points a, c and e are impact nodes whereas b and d are cross 

# 

nodes. Let (Z ,Z ) be the state of the loose mass. Here onwards, 

P P 

the subscripts + and — are used, respectively, to denote the time 
instants just after and just before the crossing of a node. The 
motion of the system is given by 
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Figure 5.3 Half cycle of a typical four impacts per cycle motion. 
Arrows indicate the locations of impact. 
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i+i 

("C) = exp(-hj_^^ (T-x^) ) (A^ cos t)(x-t^) + sin 7){x-r^)) 

(5.4a) 

and 

• 1+1 

Z (x) = constant for ( x ) , < x < ( x . ) 

^ i 

(5.4b) 


where V - 

denotes the instant corresponding to the i^ and twice of h^^^ 
is the value of the function c.f between the i^^ and (i+1)^^ 
nodes . 

If the i^^ node is an impact node, then 



♦ • • 

Zl(\)- = “i Zp(^)- “2 + 

♦ ♦ • 

= “3 •" “4 

where 

ocj^ = (X - ) / ( 1 + X) 

oc^ = ( 1 + ) / ( 1 + X ) 

= X ( 1 + ) / (!+:»:) 

= (1 - X ) / ( 1 + X ) • 

Further, at an impact node, 


The symmetry of the motion implies 

and since the velocity of the loose mass remains unchanged between 
two successive collisions, we get 
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2p(S)- = - (5.9) 

for the four impacts per cycle motion as depicted in Figure 5.3. 

+* Vi 

If the i node is a cross-node, then the following conditions 
hold good : 


Zl('Ci) = + 1 

and 


(5.10) 


Zl(T )- = (5.11) 

Using the above relations and lengthy algebraic manipulation, we 
get n number of equations in 's where n represents the number 
of nodes per half cycle of the motion. Without loosing any 
generality, we can set x^ = 0 and the remaining (n-1) 

transcendental equations can be numerically solved for x^ ^s. 
Thereafter the other unknowns, viz., A^'s, s and the velocity 

of the loose mass between two consecutive impacts, can easily be 
obtained. 

The stability of the solution obtained above can be analysed 
by the method of error propagation [68]. A suitably chosen initial 
disturbance vector Ci having four components is introduced into 
the analytically constructed solution . Let C.j. the disturbance 
vector measured after a cycle , then related to C.j. 


= [P] + higher order terms . (5.12) 
The matrix [P] can be estimated in closed form, for a piece-wise 
linear system, as described below. 

Let the i^^ node be an impact node. In that situation, as 
revealed by the numerical solution, the (i+1)^^ node is a cross 
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till 

node and the (i+2) node is again an impact node and so on. 
til til 

Between the i and (i+1) nodes, the motion of the system is 
given by .equations (5.4). A lengthy algebraic manipulation finally 
gives 

= Q. Zi(\) + h Zi(\)+ . 

* 

and = G Z^(T^) + U Z^(r^)^ 

where E , D , Q , L , G and U are functions of t , r . , and 

i' i' *i' l' 1 1 i' l+l 

r etc. Since in between the impacts the loose mass moves with a 
1+2 

constant velocity, we have 


(5.13) 

(5.14) 


w + = 


^ 1+2 ” 


(5.16) 


Introducing a perturbation in the i node given by the vector 
Cj = ( At , AZ , AZ^ , AZp ) ^ and following this perturbation 
through equations (5.4)— (5.16) , the following linearized relation 
is obtained ; 


^ 1+2 




(5.17) 


Where is a 4X4 matrix. Thus, finally we can write 
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1+4 


= cp]!I! c 


1+2 ''i+2 ’• ■* 1+2 


Cj = [p] c. 


(5.18) 


So far we have considered only one-half of the cycle, but the 
motion being symmetric, the other half behaves in a similar 
fashion and the matrix [P] can be obtained analytically. 

If the eigenvalues of [P] lie within the unit circle about 
the origin of the complex plane, then the corresponding solution 
is stable. When a control parameter is varied, the eigenvalues may 
cross the unit circle in different ways leading to different types 
of bifurcations of co-dimension one. 


5.4 Numerical Method of Stability Analysis 
The analytical method presented in section 5.3 becomes 
unwieldy with increasing number of collisions per cycle. 
Therefore, a numerical method is necessary for obtaining both 
the solution and its stability analysis. A solution obtained by 
numerical integration is of course stable. However, nothing is 
known about the local invariant subspaces (eigen structure) of the 
solution. In this section, a method for calculating the eigen 
properties of the fixed point of the Poincare map of a numerically 
obtained solution is presented. 

Let a Poincare section be defined as 

Sp = { (T, Z^, 2 ^, 22 ) : 2^ = r, \> 0 } (5.19) 

where y is a constant. A fixed point W in corresponds to a 
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periodic solution and is defined as 


* * 

W = P ( W ) 


(5.20) 


Neither P, nor W can be calculated analytically. For ascertaining 
the local nature of the fixed point, 


DP = 


[ 4 ] 


(5.21) 


with 


W = { T , Z^, Z^} 


needs to be computed. As the system under consideration is 
autonomous, DP can be written in the following partitioned form: 

1 


DP = 


0 

0 

0 


^1 ^2 ^3 


[B] 


(5.22) 


Where [B] is a 3x3 matrix and is defined as 

[B] = [ b^j ] ; i,j = 1 ,2 and 3 (5.23) 
Let us assume that P can be expanded in a Taylor series around the 
fixed point W*. In what follows, we explain how [B] can be 
constructed. Introducing a series of perturbations. 




' 1 ' 

AZ., 

- = A ^ 

0 > 

2 

AZ- 

V 2 J 

s 

0 


(5.24) 


where 


s 


. . . N , with N representing the number of 
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perturbations used, in the numerically obtained (integrating 
equations (5.1)) fixed point, their new values are obtained after 
the first return on as {AZ^ ,Az/ , Az/}^ . One can write 

‘ / N-1 n 

AZ. = Ag ^ ^n ^ ^s ^ where s = 1,. . . N (5.25) 

n=l 

Thus, N simultaneous equations in , and X^'s are obtained which 

can be solved to get One can go on increasing the value of N 

to get a converged value of b^^. However, if A^ is taken 

sufficiently small, a very low value (around 3 or 4) of N is 

needed for convergence. However, too small a value of A is also 

s 

discouraged in order to avoid computational errors that occur 

while handling small numbers. Other elements of [B] can also be 

similarly calculated. For the elements of the first column, we 

need perturbations like A { 1 0 0 }^. Similarly, for the second 

and third columns we use A { 0 1 0 and A { 0 0 1 }^, 

s s 

respectively. Once DP is constructed numerically, then the 
movements of the eigenvalues of DP near the boundary of the unit 
circle give an indication of the nature of the impending 
bifurcation. 


5.5 Results and Discussions 

It may be noted that the parameter d^ characterizes the 
non-linearity of the damper. Thus, it is pertinent to analyse the 
dynamic transitions under the control of the parameter d^. The 
whole discussion has been divided into two parts. In the first 
part, the existence , stability and bifurcations of the basic 




periodic solutions with higher (> 2) number of impacts are 
discussed. The nature of the interwindow chaos occurring between 
two basic windows are also analysed. Unlike the motions with 
higher number of impacts, the solutions consisting of two impacts 
per cycle show several distinct bifurcational behaviours which are 
very sensitive to the other parameter values. Therefore, a 
separate treatment of this type of motions is required. The second 
part is directed to a detailed discussion of the stability and 
bifurcations of the periodic solutions with such two symmetric 
impacts per cycle. 

5.5.1 Bifurcations of the Periodic Solutions with More than Two 
Impacts 

In this section, the stability of the basic periodic solution 
with four impacts per cycle for the piecewise constant 
self-exciting force is investigated. The periodic solutions with 
higher number of impacts and other continuous self-exciting 
functions are treated numerically. 

5.5. 1.1 Analytical Results 

The bifurcation structure of the symmetric motion with four 
impacts per cycle is studied in the (x, parameter space. The 
stability boundaries in the parameter space are delineated in 
Figure 5.4. The eigenvalues, at both the right and left stability 
boundaries, are seen to cross the unit circle along the positive 
real axis. Beyond the right stability boundary, the non-existence 
of any solution implies the occurrence of a saddle-node 


Figure 5.4 Stability boundaries of symmetric four impacts per 

cycle solution. - - - saddle-node Jboundary, 

pitch-fork boundary, c = 0.07, r = 0.1. 
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jjj^furcation . This saddle-node bifurcation is best depicted in 
Figure 5.5, where the limit cycle amplitude is plotted against the 
control parameter d^. The details of the consequence of such a 
bifurcation are discussed later. An unstable solution continuing 
beyond the left stability boundary (shown in Figure 5.4) implies 
the occurrence of a pitchfork (symmetry-breaking) bifurcation. The 
consequence of a pitchfork bifurcation (supercritical) is the 
evolution of an asymmetric solution with four impacts per cycle. 

5.5. 1.2 Numerical Results 

For numerically integrating the equation of motion, a 
fourth order Runge-Kutta-Merson algorithm (NAG routine) with 
variable step size control has been used. This numerical scheme 
has also been used for zero contact-time model (Figure 5.1(a)) 
which has been treated analytically in section 5.3.1. The 
amplitudes of the limit cycle obtained (a^) analytically and 
numerically are shown in Figure 5.5. It may be concluded from this 
figure that equations (5.1a) and (5.1b) with a high value of r^(= 
20000 in- the present study) can adequately represent the impact 
phenomenon. All subsequent results are obtained by numerically 
integrating equations (5.1) which allow for a finite contact time 
during the collision. As mentioned earlier, the numerical 
simulation revealed the existence of several windows of periodic 
solution, each with a specific number of collisions per cycle, in 
the X ~ plane. As shown in Figure 5.6, these windows are 
separated by zones where the oscillation undergoes a series of 


bifurcations. 



Figure 5.5 Bifurcation diagram of symmetric four impacts per 

cycle motion. =0.39. e = 0.07. r^= 0.1. stable 

solution. unstable solution. SN - saddle-node 

instability, PF - pitch-fork instability, o - results 
obtained via numerically integrating equation 5.1 for 

rj= 20000. 




Figure 5.6 Stability boundaries for r^= 0.1 with the form of f as 
shown in Figure a) 5.2(d), e = 0.07; b) 5.2(b), e = 
0.1; c) 5.2(c), e = 0.15; d) 5.2(a), e = 0.15. w6 
denotes the window of six impacts per cycle symmetric 
motion, w4 denotes the window of four impacts per 
cycle symmetric motion. - - - flip bifurcation line. • 
symmetric period-three solution; (•) asymmetric 
period-three solution. 
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The periodic solution of sa.ch window undergoes a transition 
to chaos via intermittency when d^ is increased beyond a critical 
value (for a fixed value of x) • The intermittency has been 
identified by its typical signature of a long, laminar, almost 
regular phase of oscillation interrupted by relatively short 
chaotic bursts as shown in Figure 5.7. 

A nice phenomenological description of the origin of 
intermittency can be found in the references [115-117]. This 
description was with particular reference to non-linear maps which 
are equivalent to a Poincare map of the dynamic evolution of a 
system described by differential equations. Three different types 
of intermittency, based on the local bifurcation of a periodic 
solution, have so far been identified. Mathematically, type I 
intermittency is associated with a saddle-node (tangent) 
bifurcation, whereas types II and III are preceded, respectively, 
by subcritical Hopf and flip bifurcations. 

In the analytical study presented in section 5.5.1, a 
saddle-node bifurcation was revealed at the right boundary of the 
window of the four impacts per cycle symmetric solution. The 
occurrence of a saddle-node bifurcation at the right boundary of 
every wihdow, independent of the form of function f, has been 
confirmed by the numerical method outlined in section 5.4. So the 
ensuing intermittency is classified as type I. 

Another way of testing a type I intermittency is to plot a 
Poincare map and detect the well-known narrow channeling effect 
[117]. A Poincare map shown in Figure 5.8 clearly demonstrates 
this narrow channeling effect. 






142 


Another set of Poincare map tt o Z , with n denoting the 

# 

projection on ( '^>2 ) plane, has been plotted in Figure 5.9. 
These plots clearly show a high measure density (in some measure 
theoretic sense) near the ghost of the fixed point which existed 
prior to the bifurcation and a low measure density elsewhere. The 
points with a high measure density are associated with the local 
dynamics (laminar phase) whereas those with a low measure density 
are related to the global excursion i.e., the chaotic bursts. It 
should be mentioned that the structure of the attractor shown in 
Figure 5.9 is insensitive to the models of the self-exciting 
mechanism. Somewhat detailed study of this chaotic attractor will 
be presented later . 

The statistical properties associated with the intermittent 
signature have also been calculated. The following scaling law is 
obtained : 


<^> - c“ ^ , 5 ~ 0.52 
c 

c c » 

where <l> is the average laminar length and = d^ - d^ and d^ is 
the critical value of the parameter d^ at the onset of 
intermittency. The above scaling law is very close to what has 
been predicted analytically (5 = 0.5) by Pomeau et al.[117]. If d^ 
is increased beyond the right boundary of a window (see Figure 
5.6), the intermittent oscillation encounters frequent chaotic 
bursts giving rise to chaos. Some chaotic attractors are shown in 
Figure 5.10. The shapes of these attractors are very similar to 
those shown in Figures 5.9. The only difference is that the points 




Figure 5.9 Poincare sections corresponding to intermittency. x - 

0.45, r = 0.1, r = 0.5. a) f - as in Figure 5.2(d), e 
m 

* 0.07, d = 0.736; b) f - as in Figure 5.2(b), e = 
o 

0.1, d = 0.773; c) f- as in Figure 5.2(c), e * 0.15 , 
o 

d^ss 0.532; d) f- as in Figure 5.2(a), e ■ 0.15, d^» 
0.4505. 





Figure 5.10 Poincare sections corresponding to chaos ensuing 
intermittency. =* 0.45, r^ = 0,1, y = 0.5. a) f - as 
in Figure 5.2(d), e * 0.07, d^ = 0.76 ; b) f - as in 
Figure S.2(b), e » 0.1, » 0.79; c) f- as in Figure 

5.2(c). c * 0.15 , d^« 0.55; d) f- as in Figure 
5.2(a). c » 0.15, d - 0.46. 
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on the attractors in Figure 5.10 are more evenly distributed. 
Thus the saddle-node bifurcation followed by type I intermittency 
is seen to be persistent and is independent of the details of the 
form of the function f. 

It has already been noticed in section 5. 5. 1.1, that the left 
stability boundary of the symmetric solution with four impacts per 
cycle corresponds to a pitchfork bifurcation. The numerical 
solution has revealed that this feature is also true for all the 
windows and for all variations of the function f. Two asymmetric 
solutions follow this pitchfork bifurcation. These asymmetric 
solutions have the same number of impacts per cycle as that within 
the window under consideration. With a further decrease in d^, the 
asymmetric solutions undergo a flip bifurcation giving rise to a 
period-doubling sequence (see Figure 5.6). A further decrease in 
leads the period-doubling sequence into chaos. The whole 
dynamical transition is shown in Figure 5.11. 

It has already been mentioned that there exist a whole lot of 
complicated solutions between two windows of basic solutions. One 
of these solutions is a period-three oscillation. Just like the 
window of symmetric basic solutions, a window of symmetric, 
period-three solution is found to exist irrespective of the form 
of the function f. This window is very narrow and therefore, only 
a set of representative points indicating the location of the 
window is shown in Figure 5.6. These symmetric period-three 
oscillations undergo bifurcations similar to those shown by the 
basic solutions. Type I intermittency is observed at the right 
boundary of the window as shown in Figure 5.12. At the left 



Figure 5. II Period doubling bifurcation to chaos in poincare 

plane with ^ = 0.5. f - as shown in Figure 5.2(a), e 

= 0.15 ,r = 0.1 , = 0.45 . a) d = 0.36, b) d = 

m o o 

0.35 , c) d = 0.338 , d) d = 0.334, o and A denote 
o o 

two different solutions. 







ZjCi) 


Figure 5.12 Intermittency of symmetric period-three oscillation 
in Poincare plane with r = 0.5. f- as shown in Figure 
5.2(a), e = 0.15, x = 0.42, - 0.1 , = 0.3211. o 

- symmetric period-three oscillation prior to 
intermittency at d_ ■ 0.321 
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boundary, symmetry-breaking takes place giving birth to a pair of 
asymmetric period-three oscillations which ultimately lead to 
chaos via period-doubling bifurcations. A transition of this kind 
is shown in Figure 5.13. 

Besides the symmetric period-three window, an asymmetric 
period-three window is also seen to exist whose location is again 
identified by a set of points in Figure 5.6. With decreasing 
values of d^ for a fixed x, such a pair of asymmetric solutions 
undergo period-doubling bifurcations and ultimately lead to chaos. 
This transition is shown in Figure 5.14. 

The existence of an intermittency route to chaos immediately 
after a saddle-node bifurcation was mentioned in the beginning of 
this section. To analyse the structure of the attractor 
corresponding to this intermittency, let us consider Figures 
5.15(a) and 5.15(b). The existence of the region L in Figure 
5.15(a) with a high measure density of the return points can be 
explained by the existence of a saddle-node prior to the 
bifurcation. But the region G with a low measure density of the 
return points corresponds to some global phenomena. It may be 
concluded that this region is the outcome of subharmonic 
instability, either of the basic solution (see Figure 5.11) or of 
the period-three solution (see Figures 5.13 and 5.14). Mostly it 
was found that a period- three solution exists near the saddle-node 
boundary and the region with a low measure density is related to 
the subharmonic instability of this solution. 
















5.5.2 Bifurcations of the Periodic Solutions with Two Symmetric 
Impacts 

So far we have discussed the existence and stability of 
the periodic solutions with more than two impacts. It is seen that 
for a fixed combination of {Xt and c ) , as d^ is increased, 
several windows of periodic solutions (having the basic period) 
arise. Each of these windows corresponds to a specific number of 
impacts per cycle . In fact, the number of impacts per cycle 
corresponding to this basic period is seen to decrease with 
increasing d^. Ultimately, if the other parameters are suitably 
chosen, a symmetric motion with two impacts per cycle appears. In 
real life applications of an impact damper, this mode of operation 
is most effective. As suggested by the preliminary numerical 
simulation, a periodic motion with two impacts per cycle is seen 
to bifurcate in a number of ways depending on the values of r^^^, x 
and e. In what follows, we analyse the bifurcations of this 
solution and try to show the qualitative equivalence, if any, for 
all the forms of self-excitation considered hitherto. 

5.5.2. 1 Analytical Results 

To bring out the essential dynamical features of the 

system, the bifurcational behaviour are studied under the control 

of the parameter d^ with fixed values of e and x and several 

choices of the value of r . Different bifurcational sequences are 

HI 

observed depending on the value of r^^^. The bifurcation diagrams 
are plotted in the (a^ vs d^) plane as shown in Figure 5.16. 
Following observations can be made from this figure: 



Figure 5.16 Bifurcation diagram of the periodic motion with two 

symmetric impacts per cycle for various forms of f a) 

as shown in figure 5.2(d); e =0.2, x =0.2, for 

r = 0.1, for r = 0.12; b) as shown in figure 

m m 

5.2(b); X =0.2, e = 0.25, for r^ = 0.1, 

for r = 0.14; c) as shown in figure 5.2(c); c = 0.3, 
m 

X - 0.2, for r = 0.12, for r = 0.14; d) 

as shown in figure 5.2(a); c = 0.3, x = 0.2, 

for r = 0.1, - - - for r = 0.15. PF- 

m .m 

pitchfork, SN- saddle-node, HB- Hopf bifurcation. 
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(i) For fixed values of e (<< 1) and x {< 1 ) and a range of values 
of r^, a window of periodic solutions with two impacts per cycle 
exists. 

(ii) For low values of r^, as d^ increases beyond a critical 
value, a symmetric solution with two impacts per cycle is 
generated via a supercritical pitchfork bifurcation of its 
asymmetric counterparts. This symmetric solution persists for a 
range of values of d^. Beyond another critical value of d^, the 
symmetric solution loses its stability via a supercritical Hopf 
bifurcation. 

(iii) For somewhat higher values of r^^^, multiple solutions are 
resulted. The solution branch of higher amplitude is generated by 
a pitchfork bifurcation and terminated by a saddle-node 
bifurcation. The solution branch with lower amplitude undergoes a 
saddle-node bifurcation at the lower end of d^ and a supercritical 
Hopf bifurcation at the higher end of d^. 

(iv) For still higher values of r^^, the solutions with two impacts 
per cycle become either unstable or non-existent. 

Cv) The bifurcational features of the periodic solutions are seen 
to be qualitatively equivalent for all types of self-exciting 
function considered in this paper. 

It may be noted that the above mentioned scenario of the 
bifurcational behaviour can also be seen by varying c with fixed 
values of r and x- 
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5. 5 . 2. 2 Numerical Results 

The local stability analysis, discussed in the previous 
section, has revealed different kinds of bifurcations of a 
symmetric periodic solution with two impacts per cycle namely, 
saddle-node , Hopf, pitchfork etc. In this section, numerical 
results are presented to first unfold the dynamic evolution of the 
system beyond the above mentioned Hopf bifurcations. These results 
are shown in terms of Poincare sections in Figures 5.17-5.19 for 
the three different types of self-exciting functions shown in 
Figures 5.2(a)-5.2(c) , respectively. 

The Hopf bifurcation of a symmetric solution with two impacts 
per cycle gives rise to a quasiperiodic attractor which is 
identified by a closed curve in the Poincare section Z . It is to 
be noted that the quasiperiodic attractor is smooth in nature near 
the bifurcation point (Figures 5 . 17 (a) , 5. 18 (a) , 5 . 19 (a) ) . However, 
when the control parameter d^ is increased further, the smoothness 
of the attractor is destroyed (Figures 5.17(b), 5.18(b), 
5.19(b),). Such non-smooth attractors are typical of dynamical 
systems having piecewise smooth characteristics and have been 
studied by several investigators [41]. With further increase in 
dg, the quasiperiodic motion gets locked into a periodic attractor 
of higher (than the basic) period which subsequently becomes 
unstable and chaotic. Such a quasiperiodic route to chaos via 
frequency locking is apparent in Figures 5.17-5.19. This 
particular route to chaos is found to be qualitatively very much 
in accordance with the dynamics of Curry-Yorke model and several 










Figxire 5.18 Quasiperiodic route to chaos in Poincare plane for 
the form of the function shown in Figure 5.2(b) with 
various values of d^. y =0.1, e =0.15, r^ = 0.075, x 
= 0.2. a) d„ = 2.5, b) d^ = 2.6, c) d = 2.68, d) d^ 









Quasiperiodic route to chaos in Poincare plane for 
the form of the function shown in Figure 5.2(c) with 
various values of d^. T ~ O-)* “ 0.14, x ~ 

0.2. a) d„ = 1.5, b) d- = 1.6, c) d = 1.653, d) d^ = 


Figure 5.19 
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other experimental observations of hydrodynamic instabilities 
[117]. 

Now let us return to the numerical results which are 
analogous to those shown in Figure 5.15. Here | instead of 
a^, is plotted again; (Figures 5. 20 (a) -5 . 20 (d) ) for different 
values of r . At low values of r , only one value of j Z. j is 
observed (Figure 5.20(a)). However, for somewhat higher values of 
r^^, multiple solutions are obtained. Depending on the location of 
the saddle-node points, hysteresis involves jump between different 
types of solutions as depicted in these figures. For still higher 
values of r^, the periodic solutions with a low amplitude may be 
completely destabilized and are replaced by quasiperiodic 
solutions or other solutions bifurcating from the same (Figure 
5.20(b)). If the value of r^ is increased further, then the lower 
amplitude solution ceases to exist while the other periodic 
solution with two impacts per cycle continues. This solution, 
beyond a critical value of d^, is terminated by a saddle-node 
bifurcation giving rise to type I intermittency (Figure 5.20(c)). 
Ultimately, at some high value of r^, periodic solutions with two 
impacts per cycle cease to exist (Figure 5.20(d)). 

The delineation of regions of periodic solutions with four 
impacts per cycle and two impacts per cycle are apparent in 
Figures 5.20(a)-5.20(d) . Complicated transitions take place in the 
intervening region. We have already discussed similar transitions 
with reference to the interwindow space between six and four 
impacts per cycle. 
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Figure 5.20 Bifurcation diagram of the periodic solutions with 

two symmetric impacts per cycle for the form of the 

function f shown in Figure 5.2(b) with various values 

of r . y = 0.4, e * 0.25. a) o two impacts per cycle 
m 

and • four impacts per cycle for r^= 0.05; A two 
impacts per cycle and A four impacts per cycle for r^ 
= 0.075; o quasiperiodic motion; b) o two impacts per 
cycle and • four impcts per cycle for r^ = 0.08; ® 
two impacts per cycle for r^ = 0.085; o quasiperiodic 
motion; c) o two impcts per cycle and • four impacts 
per cycle for r^ = 0.09; d) • four impacts per cycle 
, no symmetric two impacts per cycle exists for r^ = 
0.18. IW represents interwindow space. 





5.6 Conclusions 


Dynamics and. bifurcations of a class of autonomous 
self-excited systems with momentum transfer mechanism have been 
studied in detail. Both analytical and numerical results are 
presented. The chosen parameter space is shown to be spanned by 
several windows of periodic solutions of a basic period with the 
interwindow spaces filled up by complicated motions. Each window 
of periodic solutions corresponds to a specific number of impacts 
per cycle. The number of impacts per cycle corresponding to a 
window of basic solution are seen to decrease with increasing 
dp. The periodic solutions with higher number of impacts are seen 
to undergo type I intermittency transition to chaos following a 
saddle-node bifurcation when a control parameter is increased 
beyond a critical value. However, the periodic solutions with two 
symmetric impacts per cycle undergo jump or intermittency 
(depending on the values of other parameters) after a saddle-node 
bifurcation . A quasiperiodic route to chaos is preceded by a Hopf 
bifurcation associated with a symmetric motion with two impacts 
per cycle. When the control parameter is decreased below a 
critical value, each symmetric solution of the basic period 
undergoes a pitchfork bifurcation giving birth to a pair of 
asymmetric solutions which ultimately generate a chaotic motion 
through a period-doubling sequence. The symmetric, period-three 
solutions undergo bifurcations which are similar to those 
exhibited by the basic solution. All the above phenomena are found 
to be independent of the model of the self-excitation mechanism. 


CHAPTER-6 


BIFURCATIONS AND CHAOS IN A HARMONICALLY 
FORCED SELF-EXCITED OSCILLATOR WITH IMPACT DAMPING 

6.1 Introduction 

A detailed account of the optimal design and performance of 
an impact damper for controlling the high-amplitude vibration of a 
forced self-excited oscillator is presented in chapter 4. It has 
been shown that an impact damper may be very effective in 
neutralising the vibration near and around the harmonic 
entrainment zone. However, no mention has been made about the 
types of instabilities that might impair the performance of such 
an impact damper. None the less, a designer should know the 
limitations of the device in its operating conditions. Moreover, 
an impact damper together with its most simplistic low dimensional 
model also represents a class of mechanical devices with loose 
sub-components which introduce clearance non-linearity and impact 
damping in the system. Owing to the presence of non-linearity such 
systems can exhibit a variety of non-linear instabilities and 
complex motions which call for an in depth study from theoretical 
as well as practical points of view. Keeping the above facts in 
mind, the different types of instabilities, bifurcations, and 
chaotic motions of a harmonically forced, single degree of freedom 
van der Pol oscillator with an impact damper are investigated. 

It has been shown in previous studies that the motion with 
two symmetric impacts per cycle is the most effective operating 
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iTiod.© of impact dsmpsir • This is tjrvis iirirsspsctivs of th© typs of 
damping ( dissipative or self-exciting) present in the primary 
oscillator* Sung et al.[l01] have studied the bifurcational 
behaviour of an impact damper for a forced, linear, strictly 
dissipative oscillator near its operating mode with two impacts 
per cycle. Several types of codimension one bifurcations viz., 
pitchfork, saddle— node and Hopf etc. have been detected. Cascades 
of period doubling bifurcations along with global bifurcations 
such as homoclinic tangle of stable and unstable manifolds have 
also been studied. Chatter jee et al. [114] have investigated the 
super structure of the bifurcational behaviour of an autonomous, 
self-excited oscillator with impact damping. However, no studies 
concerning the bifurcations of a forced self-excited oscillator 
with impact damping are known to have been reported. 


6.2 Equations of Motion and Method of Solution 

Equations of motion of a forced self-excited oscillator with 
an impact damper are already given by equations (4.1) and (4.2). 
As motions with symmetric, two impacts per cycle are relevant from 
the practical point of view, the present study is mainly centered 
around this type of motion. The EPL method described in section 
4.4.1 is used to construct the periodic solution and the method of 
error propagation elaborated in section 4.4.2 is applied to 
analyse the stability and bifurcation of the solutions. 
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6.3 Theoretical Results 

In this section we present different possible bifurcations 
obtained analytically using the EPL technique and the associated 
error"propagation method. There are several parameters in the the 
model of the system. In practice, it is very difficult to handle 
all of these simultaneously. However, as we are interested in the 
local bifurcations of codimension one, the parameters are selected 
so as to detect all the possible types of bifurcations. Several 
bifurcation sets are depicted in Figures 6.1 - 6.3 in different 
parameter spaces. The types of bifurcations (subcritical or super 
critical) are also delineated. The existence of cusp catastrophes 
are particularly apparent in Figures 6.2 and 6.3. Single (Figures 
6.2(a), 6.2(c)) and double ( Figure 6.2(b), rjjj= 0.5) cusps are 
seen in the bifurcation sets obtained in d^ - e plane. In the F^ - 
n plane (Figure 6.3), however, two (at E and F) or three (at A, B 
and C) cusps are seen to occur. A direct numerical simulation of 
the equations of motion are employed to analyse the dynamical 
bearings of these bifurcations. These results are presented in the 
following sections. 


6.4 Numerical Results 

For numerical integration, we have considered equations 4.2 
in which the impacts are modelled by springs and dashpots. A 
fourth order Runge-Kutta-Merson algorithm with adaptive step size 
control has been employed. Numerical integration revealed 
different types of intermittent and quasiperiodic routes to chaos 
as elaborated below. 


Figure 6.1 Bifurcation set in r^- plane, - I, ^ = u. /, r^- 
0.5, e = 0.1. The zone enclosed between the lines AB 
and CD corresponds to stable symmetric motion with 
two impacts per cycle. AB and CD correspond to 
supercritical and subcritical Hopf bifurcations, 
respectively. 
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Figure 6.2 Bifurcation sets in e - plane. Q - 1. a) Fq= 0.5, 

r = 0.25; b) F = 0.5, = 0.7; c) = 0.7, r = 0.25. 

m 0 

Hopf bifurcation (the upper line corresponds 

to supercritical bifurcation and the lower line 
corresponds to subcritical bifurcation), 
pitchfork bifurcation (subcritical). 
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641 Type II Intermittency 

The zone of stable solutions with symmetric two impacts per 
cycle is delineated in the r^- d^ plane as shown in Figure 6.1. 

local stability analysis reveals that both the stability 
boundaries (marked as AB and CD) correspond to Hopf bifurcation. 
AS the non-linear stability analysis turned out to be complicated, 
the form of the bifurcation could not be established analytically. 
However, it is well known that two possibilities exist when a 
periodic motion encounters a Hopf bifurcation (secondary Hopf 
bifurcation or Neimark bifurcation) . Either a quasiperiodic motion 
results (when the bifurcation is supercritical) or a more 
complicated evolution (when the bifurcation is subcritical) having 
a local-global or entirely global feature appears. In the present 
case, the boundary AB turned out to be a supercritical bifurcation 
line as the subsequent motion is strictly quasiperiodic. However, 
when a parameter is varied (here d^ is decreased while keeping r^^^ 
constant) across the boundary CD, an intermittent evolution 
appears. Such an intermittent signature is shown in Figure 6.4. A 
closer look into the signature reveals two distinct phases of the 
motion. One of these consists of frequency modulated oscillation 
(so called laminar phase) with diverging envelop. Two such phases 
are punctuated by other phases which are apparently disordered 
(chaotic phase) . According to the classical Pomeau-Manneville 
[115-117] categorization of different types of intermittencies 
based on local bifurcations, the present intermittency is of type 
II as the preceding bifurcation is a Hopf bifurcations. This 
indirectly implies that the Hopf bifurcation on the boundary CD is 
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gubcritical - The average length <l> of the laminar phases is seen 
to comply with the scaling law (Figure 6.5) predicted by 

Pomeau-Manneville and described by 

<l> ~ In (l/M) 

where h = and d^ is the bifurcation threshold. 

Here, it is worth mentioning that though the type II intermittency 
is known to exist, it rarely occurs in physical systems. Only very 
few examples are found in the literature [118-119]. 

6.4.2 Intermittency After Subcritical Pitchfork Bifurcation 

The zones of stable two impacts per cycle solution are 
delineated in Figure 6.3 in the F^ vs n plane. The boundary AB in 
this figure can be analytically shown to correspond to a 
subcritical pitchfork bifurcation. Immediately after this 
bifurcation (n is decreased across the boundary while keeping F^ 
constant) , another intermittent motion appears as shown in Figure 
6.6. In this figure, the time history of the Poincare iterates (so 
called strobe map with the strobe frequency fi) is plotted. Such a 
kind of intermittency created by a subcritical pitchfork 
bifurcation is not reported so far. 

Theoretical models of the dynamics of intermittency have been 
put forward by Pomeau-Manneville [115-117]. For describing the 
mechanism of intermittency they used non-linear maps which may 
stand for a suitable Poincare map of the periodic solutions of an 
ordinary differential equation. In their model, a kind of narrow 
channeling effect accounts for the laminar phase and a proper 
reinjection principle explains the recontinuation of the laminar 
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phase after a chaotic burst. As the whole mechanism is based on 
the availability of a suitable non-linear map, only three types of 
intermittency could be explained each corresponding to a type of 
local bifurcation viz-, saddle-node (type I), subcritical Hopf 
bifurcation (type II) and subcritical period multiplying 
bifurcation (type III) of the fixed points of the map. However, 
the pitchfork bifurcation (or symmetry-breaking) of a periodic 
solution of a differential equation does not have any counterpart 
in the bifurcation of a fixed point of a non-linear map. Perhaps 
this is why no model has been studied for the occurrence of 
intermittency following a pitchfork bifurcation. Another reason 
may be that the symmetry breaking is a non-generic bifurcation in 
the sense that, it is structurally unstable under asymmetric 
perturbations however small. Of course, for a more detailed study 
of the mechanism of this type of intermittency a model should be 
constructed and according to the above discussion such a model 
should be based on a low dimensional differential equation rather 
than a non-linear map. 

6.4.3 Quasiperiodic Route to chaos 

The boundaries corresponding to supercritical Hopf 
bifurcation are delineated in Figures 6.1 - 6.3. When a parameter 
is varied (keeping the other constant) across these boundaries, 
the symmetric, two impacts per cycle solution bifurcates to a 
quasiperiodic solution. Such quasiperiodic motions can lead to 
chaos with further variation of the control parameter beyond the 
bifurcation line. Several routes to chaos via quasiperiodicity 
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is veil known that a great deal of hydrodynamic instabilities are 
related to different kinds of intermittency catastrophes and 
transitions. The present model, a purely mechanical 
demonstrates a variety of similar instabilities. 


one, 


CHAPTER-7 

CONCLUSIONS 


7.1 Conclusions 

The major conclusions of the present thesis are listed below. 

(i) The proposed EPL technique contains information in both time 
and frequency domains and consequently gives better results than 
the HBM which is simply a frequency domain technique. This method 
is particularly suitable for obtaining the time domain response of 
a class of piecewise smooth, non-linear systems. 

(ii) An impact damper is effective, both in passive and active 
on-off forms, for controlling the vibration of non-linear systems 
having Duffing type restoring force characteristics. Both the 
primary and secondary resonances can be controlled. A simple 
formula is obtained for providing an initial design and has been 
shown how to improve upon to obtain a final optimum design. 

(iii) An impact damper can be used to control both free and forced 
self-excited vibrations. Sensitivity of the designed performance 
towards the inaccuracy of the estimated model parameters is to be 
incorporated in the design rules. An empirical design formula for 
an impact damper attached to a forced van der Pol oscillator is 
provided after a detailed parametric study. The role of the 
coefficient of restitution for various situations has been clearly 
brought out. 
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(iv) Various instabilities occur in an autonomous self-excited 
system with an impact damper. The gross bifurcation superstructure 
of the systems under consideration is more or less independent of 
the exact nature of the model of the self— excitation. A numerical 
technigue has been developed for analysing the eigen properties 
of the stable solutions. 

(v) The dynamics of a harmonically driven van der Pol oscillator 
with an impact damper is rich with different types of 
intermittency catastrophes viz., type II, type III etc. A new type 
of intermittency after a subcritical, symmetry-breaking pitchfork 
bifurcation is seen to exist in the system under consideration. 
Several instabilities of the system has a close resemblance to the 
hydrodynamic instabilities such as Rayleigh-Benard convection etc. 

7.2 Suggestions for Future Works 
On the basis of the methods developed and results obtained in 
this thesis, the following suggestions, on research to be carried 
out in future, can be made: 

(i) The suitability of the EPL method for analysing the response 
of smooth, non-linear systems should be investigated. It should be 
emphasized, that the piecewise linearization is to be carried out 
in a dynamic sense where the junction point between the equivalent 
linear segments has to be left as an unknown quantity to be 
determined. In the usual static way of piecewise linearization 
method, this junction point is assumed a priori. 
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effectiveness of an impact damper, used both in passive 
and active manners discussed in the thesis, for controlling the 
vibration _ of non-linear systems should be verified through 
guitabl® experimentation. 

(iii) 'TThe intermittency route via the subcritical, 

syiietry breaking pitchfork bifurcation reported in the thesis 

should b- e subjected to further investigations. Especially, if 
possiblajr 3 low— dimensional mathematical model for determining its 
statisticr^cal characteristics should be attempted. It may be 
jentionecSi that the standard maps, hitherto discussed in the 
literatuicrre, cannot capture such route. 
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APPENDIX-A 


The elements of the Jacobian J of the right-hand-side of equations 
( 3 . 17 ) are as follows: 


J = 
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'^21 *^22 “^23 '^24 

■^31 '^32 '^33 “^34 

■^41 “^42 '^43 "^44 


where 


11 


h 

2 


r B(b) ^ 

ri2 = cos (02 - e^) - ^ cos 8^ , 


1 . X SB(b) 


13 


20 . 


sin (e_ - ej 


1^ ab 


14 


r^B(b) 

”20 


cos (82 - 0^) 


COS0^ - r^^ B(b) cos (02-03^) 3 ^ 


21 


2 a ^0 


4 0 


22 


1 

2 a 0 


sin8^ - r^B (b) sin (02“®i^ 


■^23 “ 


r ^ cos ( 82 - 0i ) ^( b ) 

2a0 “ 3 b 


24 


r B(b) 

- Isr - '®2 - 


J S= _ _ cos 
^31 2 


(»! - S2 > -< In - ^ ' 



196 



197 


and 


^32 


ha _ 3^3 


au„ 


sin (9 _ 0 ) _ f £a_ ^ I . -as 


J_ _ = 0 , 
33 • ' 


34 2Q 


^41 = 


^ I COS ©2 + h a n sin ( 9 ^- 9^) - 

( a(j^ + I a^)cos ( 0 ^- 02 )j- 

^ |h n sin( 8 ^- 9^) - (wj + I a^) cos (9^- e^)| , 


42 2bn 


^ -|h a Q cos '*’ T (^o" ® 


l)}' 


^43 = 


- {>= - (^ + 1 ) - tj} , 


44 


II cos (e^-e^) - ( 


3a^ 

2 fib -Ibfi ^ ^® 2 “® 1 ^ 


2 bfi ®2 


f = bn^ - r B(b) (1 + i-) ha fi sin ( 03 ^- 02 ) 
^ ^ m 


- ( awj + I a^) cos (© 2 "®!) ®2 


APPENDIX-B 

The mathematical model of an impact pair is shown in Figure 
B.l The contact dynamics of the .two masses m„ and m. are 
represented by the differential equations 



’“2>'2 =2 + >^2 

(X2-X^) = 0 

(B.l) 

and 


(X^-X2) = 0 

(B.2) 


Where the dot denotes the differentiation with respect to time t. 


Substituting = X 2 - and r^^^ = ^ 2 / 10 ^ in equations (B.l) 
and (B.2) we get 
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where 
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and u = 
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2 2 


At the beginning of the contact phase we can write 
at t =0, Y^= 0, = - V 

and at the end of the contact phase, when Y^ is again zero, one 
gets 


Y = V e'^ ‘"n^/‘"d 


where ^ 2 


Thus, according to the definition of the coefficient of 
restitution (x) 
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= exp [- rr I / / i - ] 

c 

which gives ^ = — (^) 

2i/K^ [n^+ln(x)^]^^^(l+r^)^^^ 

The above model of an impact pair is seen to be satisfactory only 
for high .values of the coefficient of restitution. The model seems 
to fail for low values of x implying a high value of ^ . A high 

C 

value of damping enhances the time of contact and thus violates 
the condition of zero contact time during impact. 



APPENDIX-C 


Let us consider the system discussed in section 3.5 with x =1 
(i.e. = 0 in Figure 3.9). The time response of the primary 
oscillator, given by various methods are shown in Figure C.l. It 
is clearly seen from this figure that, the HBM while providing a 
good estimate of the maximum value of the response during a cycle 
fails to predict the time history accurately. The EPL technique, 
containing information in both time and frequency domains, is 
capable of producing the entire time history of the response very 
accurately. It is expected that this method may work very well 
even for a smooth, strongly non-linear oscillator, an aspect which 
is worth further investigations. 


201 



202 







